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PREFACE 


This  report  provides  coastal  engineers  with  documentation  that  a  wide  zone 
of  nearshore  bathymetry  responds  to  long-term  increases  in  water  level  by 
migrating  inland  with  the  receding  shoreline.  The  dimensions  of  the  zone 
affected  depend  on  the  wave  exposure.  A  simple  procedure  is  presented  for 
estimating  the  magnitude  of  shore  recession  and  the  depth  of  profile  adjust¬ 
ment  for  any  sandy  stretch  of  shore  on  the  U.S.  side  of  the  Great  lakes. 

This  report  is  based  on  a  9-year  series  of  nearshore  surveys  conducted  on 
the  eastern  shore  of  Lake  Michigan.  The  first  three  surveys  (1967,  1969,  and 
1971)  were  carried  out  by  the  U.S.  Army  Lake  Survey  as  part  of  their  shore 
processes  investigations.  The  remainder  of  the  work  was  carried  out  under 
the  sediment  hydraulic  interaction  program  of  the  U.S.  Army  Coastal  Engineer¬ 
ing  Research  Center  (CERC). 

The  report  was  prepared  by  Edward  B.  Hands,  under  the  general  supervision 
of  Dr.  C.H.  Everts,  Chief,  Engineering  Geology  Branch,  Engineering  Develop¬ 
ment  Division,  CERC.  Reviews  and  helpful  comments  from  Drs.  C.H.  Everts  and 
R.D.  Hobson  of  CERC,  and  P.  Bruun,  are  deeply  appreciated.  Dr.  W.L.  Wood  and 
J.  Pope  provided  data  used  in  the  example  problems. 

Comments  on  this  publication  are  invited. 
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approved  31  July  1945,  as  supplemented  by  Public  Law  172,  8Bth  Congress, 
approved  7  November  1963. 


Colonel,  Corps  of  Engineers 
Commander  and  Director 
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CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI)  UNITS  OF  MEASUREMENT 


U.S.  customary  units  of  measurement  used  in  this  report  can  be  converted  to 
metric  (SI)  units  as  follows: 


Multiply 

by 

To  obtain 

inches 

25.4 

millimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.852 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  x  10"3 

kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton ,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.01745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins1 

*To  obtain  Celsius  (C) 

temperature  readings 

from  Fahrenheit  (F)  readings. 

use  formula:  C  =  (5/9)  (F  -32). 

To  obtain  Kelvin  (K) 

readings,  use  formula 

:  K  =  (5/9)  (F  -32)  ♦  273.15. 

SYMBOLS  ANU  DEFINITIONS 
0  average  height  of  affected  dunes 

d  depth  of  profile  adjustment 

d  estimate  of  d  (d  =  2.1  h^) 

h^  significant  wave  height  with  a  5-year  return  period 

k  a  constant  of  proportionality  between  wave  height  and  the  depth  of 

profile  adjustment 

in  natural  logarithm 

Q  volume  sediment  flux  into  the  survey  area 

Ra  sediment  overfill  factor — the  ratio  of  sediment  volume  supplied  bv 

profile  recession  to  that  retained  after  sediment  sorting,  packing, 
and  profile  readjustment 

sg  (z)  signum  function  having  values  of:  1  for  z  >  0;  -1  for  z  <  0;  and  0 
for  z  =  0 

T  time 

Th  thickness  of  volume  change  if  spread  evenly  over  the  survey  area 

X  average  horizontal  extent  of  profile  adjustment 

x  average  horizontal  displacement  of  the  profile  and  shoreline 

Y  longshore  extent  of  survey  area 

Z  average  vertical  extent  of  adjusting  shore  profile 

z  average  change  in  elevation  of  the  water  surface 

a  the  effective  angle  of  profile  response  if  RA  =  1  and  Q  =  0  (eq.  3) 

(also  symbol  for  "is  directly  proportional  to") 

AX  profile  digitizing  interval 
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Figure  2.  Comparison  of  annual  mean  water  levels  at  ocean  and  Great 
Lakes  sites.  During  rising  lake  levels  the  shores  of  the 
Great  Lakes  may  be  submerged  more  in  a  5-  to  10-year  period 
than  most  ocean  sites  are  in  a  century.  Reversals  in  trend 
reduce  the  longer  term  effects  on  the  lakes,  while  ocean 
sites  are  exposed  to  a  slower  but  more  persistent  rate  of 
submergence. 


PREDICTION  OF  SHORE  RETREAT  AND  NEARSHORE  PROFILE 
ADJUSTMENTS  TO  RISING  WATER  LEVELS  ON  THE  GREAT  LAKES 

By 

Bdwavd  B.  Bands 
I.  INTRODUCTION 


1 .  Purpose . 

This  report  demonstrates  that  increased  shore  retreat  during  periods  of 
sustained  high  water  is  merely  the  most  visible  expression  of  a  massive 
adjustment  affecting  a  much  wider  area  offshore.  An  equilibrium  sediment 
balance  model  realistically  describes  the  ultimate,  broad  profile  response  to 
increased  water  levels.  Beach  and  nearshore  surveys  along  a  section  of  the 
eastern  shore  of  Lake  Michigan  (Fig.  1)  in  1967,  1969,  197  1,  1975,  and  1976 
provide  the  basis  for  this  analysis.  The  results  are  generalized  to  provide  a 
simple  but  rational  approach  for  estimating  the  response  of  sandy  shore  areas 
throughout  the  Great  Lakes  to  future  long-term  change  in  water  levels. 


Background. 


Because  of  variations  in  climatic  factors  within  their  regional  drainage 
basins,  the  Great  Lakes  experience  water  level  fluctuations  uncharacteristic 
of  ocean  shores  (Fig.  2).  Dry  periods  are  common  during  which  the  mean 
elevations  of  the  lakes  decline  for  many  years  in  succession.  After  these 
long  periods  of  falling  lake  levels,  it  is  easy  for  new  property  owners  to 
overestimate  the  stability  of  the  shores  and  build  structures  too  close  to  the 
lakes.  Storm  erosion  during  later  years  of  high  water  accelerates  shore 
recession  and  creates  costly  property  damage.  The  persistence  of  high  water 
conditions  for  many  years  permits  extensive  profile  adjustment  via  erosion  and 
offshore  sediment  transport. 


Recurrently  during  periods  of  extreme  shore  erosion  there  lias  been  public 
interest  in  gaining  greater  control  over  lake  level  fluctuations.  However,  a 
study  by  the  International  Great  Lakes  Level  Board  (1973)  concluded  that  regu¬ 
lation  of  the  five  Great  Lakes,  while  possible,  would  not  provide  benefits 
commensurate  with  costs;  outflows  from  Lake  Superior  and  Lake  Ontario  have 
been  controlled  since  1921  and  1958,  respectively.  Natural  variations  in  the 
water  supply  to  these  basins  are  too  large,  however,  to  maintain  constant  lake 
levels,  so  adjustments  are  made  in  the  flows  to  benefit  the  many  interests 
involved.  Reliable  estimates  of  shore  erosion  for  various  water  level  condi¬ 
tions  are  important  in  evaluating  the  impact  of  regulation  plans  on  riparian 
interests.  An  improved  understanding  is  also  needed  for  the  proper  design  of 
coastal  construction  projects  and  beach-fill  operations,  and  the  recommenda¬ 
tions  of  setback  distances,  etc. 

A  report  on  recent  changes  in  rates  of  shore  retreat  summarized  data  from 
the  same  set  of  surveys  used  here,  but  considered  only  the  changes  within  100 
meters  of  the  shoreline  (Hands,  1979).  Dates,  types  of  data  collected,  and 
reports  on  the  earlier  surveys  are  discussed  in  Section  11. 
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11.  DATA  COLLECTION 


1 .  Profile  Stations. 

Hydrographic  surveys  were  conducted  near  Pentwater  Harbor  on  the  eastern 
shore  of  Lake  Michigan  in  both  the  spring  and  fall  of  19b9  (Pig.  3).  These 
surveys  revealed  little  variation  in  nearshore  bathymetry  beyond  the  first  50 
meters  offshore.  The  formation,  migration,  and  eventual  welding  of  an  ephem¬ 
eral  coastal  bar  to  the  subaerial  beach  constituted  the  major  change  during 
these  periods  of  relatively  limited  wave  action.  However,  when  these  profiles 
were  compared  with  profiles  from  several  of  the  same  stations  2  years  earlier, 
apparent  changes  in  bathymetry  were  evident  out  to  a  depth  of  5  meters.  To 
further  investigate  thi  apparent  long-term  profile  evolution,  profile  changes 
were  monitored  in  1971,  1975,  and  197b  by  resurveying  the  10  original  stations 
(established  in  19b7  within  a  kilometer  of  the  jettied  entrance  to  Pentwater 
Harbor)  and  24  additional  stations  spread  over  an  adjacent  50  kilometers  (Fig. 

4). 

2 .  Survey  Periods  and  Earlier  Reporting. 

Profiles  were  measured  during  six  different  survey  periods  over  a  span  of 
9  years.  The  survey  periods  and  monthly  mean  lake  levels  are  shown  in  Figure 
5.  Changes  in  bathymetry  between  1 9fc>7  and  19b9  were  reported  in  Saylor  and 
Hands  (1970).  Hands  (197ba)  provided  a  description  of  the  cross-sect ional  and 
areal  geometry  of  the  longshore  bars  throughout  the  50-kilometer  reach,  as 
well  as  information  on  grain-size  variations  and  some  speculation  on  the  pro¬ 
file  adjustment  between  19b7  and  1971.  Hands  ( 1 97 b b)  compared  profile  devel¬ 
opment  through  1975  with  a  possible  relationship  between  regional  tilting  of 
the  Great  Lakes  basin  and  variations  in  historic  bluff  recession  around  the 
perimeter  of  Lake  Michigan  over  a  120-year  period.  Hands  (1979)  incorporated 
results  from  the  1976  survey  to  describe  the  effects  of  water  level  changes  on 
the  shore  and  on  the  inner  parts  of  the  profile  (±100  meters  from  the  shore¬ 
line).  The  present  report  summarizes  adjustments  of  the  wider  responding 
profile  and  recommends  procedures  for  estimating  shore  and  nearshore  changes 
likely  to  occur  in  sandy  regions  of  the  Great  Lakes  as  a  result  of  future 
variations  in  mean  water  level  elevations. 

3.  Profile  Procedures. 

In  1967  tlie  profiles  were  measured  by  winching  ashore  a  four-wheel  level¬ 
ing  cart,  halted  every  5  meters  so  that  elevations  could  be  determined,  using 
an  engineer's  level  located  on  the  shore.  Upon  reaching  the  shore,  the  cart 
was  pulled  by  Jeep  down  the  beach  to  the  next  station  and  then  towed  offshore 
by  boat.  This  method  limited  coverage  to  depths  less  than  5  meters  and 
required  a  moderately  wide,  unobstructed  beach  for  efficient  operations.  In 
subsequent  years,  echo  sounding  was  used  to  extend  coverage  on  the  outer  part 
of  the  profile  but  instrument  leveling  continued  to  be  used  in  shallow  water 
to  provide  an  overlap  with  the  sounding  record  and  extend  the  profile  into  the 
dunes.  Boat  positioning  was  accomp L islied  by  an  optical  intersection  using  two 
transits  in  1969  and  1971  (Hands,  197ba)  and  by  a  range  and  azimuth  microwave 
system  in  1975.  In  all  years  a  transit  was  locked  on  the  profile  azimuth  for 
the  individual  station  being  sounded;  radio  contact  between  a  transit  man  and 
the  boat  operator  ensured  that  the  boat  remained  on  line. 
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periods  on  Lake  Michigan. 


All  shore  profiling  was  done  with  an  automatic  engineer's  level.  Dis¬ 
tances  were  determined  in  previous  years  by  a  rodman  carrying  one  end  of  a 

marked  measuring  wire  which  was  spooled  out  and  read  from  the  reference 

monument;  the  standard  three-wire  reading  method  was  used  in  1976.  Shore 

monuments  at  each  profile  station  were  tied  to  one  another,  to  surrounding 
bench  marks,  and  to  second-order  geodetic  monuments  surveyed  by  National  Ocean 
Survey  (NOS)  in  1973.  Vertical  reference  was  supplemented  during  profiling 

operations  using  a  system  of  water  level  recorders,  water  surface  rod- 
readings,  and  a  portable  stilling  well  which  was  placed  near  the  shoreline  at 
the  station  being  sounded. 

As  mentioned,  use  of  the  leveling  cart  limited  coverage  to  depths  of  less 
than  5  meters  in  1967.  The  outer  limit  was  extended  to  11,  16,  and  21  meters 
in  1969,  1971,  and  1975,  respectively.  No  echo  sounding  was  done  in  1 97 1> ;  the 
shore  profiles  terminated  in  about  1.5  meters  of  water. 

Ill.  PRO TILE  CHANGES 

1 .  Shore  Retreat. 

The  annual  mean  surface  elevations  of  Lake  Michigan  rose  1.4  meters  from  a 
record  low  in  1964  to  a  record  peak  for  this  century  in  1973.  The  earliest 
shore  profiles  in  the  study  area  were  surveyed  in  1967  after  the  water  level 
rise  was  well  underway.  The  rates  of  shore  retreat  from  1967  through  the  peak 

water  year,  and  for  3  years  thereafter,  are  contrasted  with  historic  retreat 

rates  by  Hands  (1979).  The  average  rate  of  shore  retreat  (landward  displace¬ 
ment  of  the  Stillwater  level)  during  the  latter  part  of  the  recent  period  of 
rising  water  was  about  six  times  greater  than  it  had  been  during  the  preceding 
120-year  period,  or  about  eight  times  greater  than  during  the  previous  50 

years.  This  increase  reflects  the  effect  of  recent  high  lake  levels.  As  the 
lake  levels  rose  the  shore  retreated  roughly  in  proportion  to  the  increase  in 
lake  levels.  Retreat  rates  remained  high  for  several  years  after  lake  levels 
stabilized;  then  as  levels  declined  between  1975  and  1976  the  beach  began 

prograding  lakeward.  During  the  last  year  of  study,  the  average  advance  of 
the  shore  was  similarly  proportioned  to  the  drop  in  lake  level  during  that 
period.  The  horizontal  change  in  shore  positions  averaged  about  4(J  times  the 
vertical  change  in  water  level  surface  during  those  same  periods.  Simple 
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linear  regression  of  shore  retreat  against  the  change  in  lake  level  explained 
50  percent  of  the  variance  in  retreat  Measurements. 

By  1975  the  shore  had  retreated  an  average  of  24  meters  from  its  1  t>  7 
position,  but  variations  between  adjacent  stations  were  large.  The  maximum 
difference  was  observed  at  Little  Sable  Point  which  lost  36  meters  in  6  years; 
the  loss  was  only  6  meters  just  2  kilometers  away.  Shore  losses  in  the 
vicinity  of  the  Pentwater  jetties  were  generally  low  due  to  a  combination  of 
shore  protection  practices  in  that  area.  Variations  among  the  other  stations 
were  not  as  easily  explained. 

More  than  80  percent  of  the  ultimate  shore  retreat  was  due  to  actual  reces¬ 
sion  caused  by  erosion  and  less  than  20  percent  was  due  to  the  immediate  effect 
of  encroachment  of  the  high  water  across  the  sloping  beach. 

2.  A  Qualitative  Description  of  Nearshore  Adjustment. 

Assume  that  an  increase  in  water  level  sets  the  stage  for  an  adjustment  of 
the  shore  profile.  The  profile  will  tend  to  follow  the  rising  water  level  by 
moving  upward  and  landward  as  the  shore  retreats.  The  zone  affected  will 
extend  from  the  point  of  highest  wave  attack  down  to  some  point  of  profile  clo¬ 
sure,  below  which  the  bottom  is  not  actively  shaped  by  surface-related  forces. 
The  point  of  profile  closure  may  be  close  to  shore  if  the  profile  is  responding 
to  a  diurnal  change  in  water  level.  However,  if  the  increase  in  water  level 
persists  for  several  years,  then  occurrence  of  the  normal  series  of  storms  may 
extend  the  point  of  profile  closure  to  depths  of  more  than  10  meters. 

Along  almost  the  entire  eastern  shore  of  Lake  Michigan,  and  at  many  other 
sites  on  the  Great  Lakes  where  there  is  sufficient  sand,  littoral  forces  have 
built  a  sequence  of  submerged  sand  ridges  or  longshore  bars  from  shallow 
inshore  to  deeper  offshore  (Fig.  6).  In  the  present  area  of  study  the  multiple 
bar  formation  extends  from  shore  to  a  depth  of  about  8  meters.  Thus,  many 
aspects  of  the  long-term  profile  adjustment  can  be  described  in  terms  of 

changing  bar  positions. 

3.  Bar  Geometry. 

Bars  in  the  Great  Lakes  have  greater  longshore  continuity  and  are  more 

regular  in  cross  section  than  those  on  most  ocean  coasts.  On  the  lakes,  long¬ 
shore  bars  are  also  persistent  from  year  to  year,  whereas  they  may  occur  only 
seasonally  on  ocean  beaches.  The  continuity,  regularity,  and  persistence  of 
longshore  bars  are  likewise  remarkable  on  enclosed  seas  (e.g.,  the  Baltic,  see 
Hartnack,  L924;  the  Mediterranean,  see  King  and  Williams,  1949;  the  Caspian, 

see  Knaps,  1966).  These  differences  probably  reflect  the  restricted  range  of 
wave  conditions  (period,  direction,  and  height)  and  tidal  variations  on  the 

lakes  and  enclosed  seas. 

Four  to  five  bars  are  persistent  from  year  to  year  at  most  stations  in  the 
study  area.  An  additional  smaller  emphemeral  sand  ridge  often  forms  closer  to 
the  shore  during  higher  wave  action,  but  migrates  to  shore  and  merges  with  the 
upper  beach  face  as  wave  conditions  wane.  In  the  longshore  direction,  these 
ephemeral  coastal  bars  are  less  continuous  than  the  outer  longshore  bars.  The 
coastal  bar  can  be  short  (less  than  1  kilometer)  and  discontinuous,  or  shore- 
tied  at  both  ends,  irregular  or  part  of  a  cellular  pattern  in  the  nearshore 
bathymetry  (Hands,  1976a).  Where  the  coastal  bar  ties  to  the  shore  there  is 
usually  a  protrusion  of  the  shoreline  and  a  flanking  indentation  (Hands,  1979). 
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Figure  6.  Longshore  continuity  of  bars  in  the  Oreat  Lakes.  The  inner  two  to 
three  bars  are  usually  visible  when  viewed  from  high  bluffs  and 

dunes  along  the  shore.  In  the  above  photo  (taken  about  900  meters 

above  lake  level)  the  inner  three  bars  can  be  seen  following  the 
curve  of  the  shoreline  from  station  10  toward  1’entwater  Harbor  in 
the  upper  right  corner  (a  distance  of  about  t»  kilometers). 

The  continuity  of  the  longshore  bars  is  interrupted  in  the  northern  part 
of  the  study  area  by  the  He ut water  Harbor  jetties.  Lai  h  year  the  outer  bars 
extend  into  the  entrance  channel  beyond  the  end  ol  the  jetties.  Typically, 
40,000  to  60,001.)  cubic  meters  of  sand  is  dredged  annually  I  rom  this  entrant  e 
bay  and  from  the  Inland  channel  where  windblown  sand  makes  an  important  con¬ 
tribution  (Seelig  and  Sorensen,  1976).  In  recent  years  an  increasing  annum 
of  ttiis  sand  lias  been  used  to  nourish  adjacent  beaches;  However,  must  cl  u  is 
taken  about  1.3  kilometers  offshore  and  dumped  into  1L  to  l 4-net et  depths. 

In  general,  jetties  which  penetrate  the  sort  zone  interrupt  the  u..i out  i 
longshore  transport  ot  littoral  dritt.  There  is  treijuent  uiiuern  turn 
jetties  divert  some  of  tile  drift  offshore  wiiere  it  ai  ■  cornu  1  a  tes  it.  water  so 

deep  that  the  sand  is  essentially  lost  t  rum  the  littoral  system.  lhv  hroad 
mour.d  of  sediment  opposite  the  Pentwater  jetties  lies  at  a  depth  of  about  i  i 

meters  and  may  have  originated  as  a  result  of  such  a  diversion  ot  longshore 

currents.  However,  the  broadly  symmetrical  appearance  and  position  ot  tuc 
mound  (Fig.  7)  suggest  it  is  more  likely  an  expression  ot  the  open  water 

uisposal  of  the  material  dredged  from  the  channel. 

An  interruption  of  longshore  bar  continuity  also  occurs  opposite  Little 
Sable  Point.  The  bars  at  this  location  are  not  only  discontinuous,  as  ob¬ 

served  in  aerial  photography,  but  are  also  much  less  regular  and  less  smooth 
in  cross  section.  Bathymetrical  ly,  Little  Sable  Point  is  a  truisitioii  /one 

dividing  the  study  area  into  two  nearly  equal  stretches  with  distimlK 
different  bar  geometry.  A  sequence  ot  tour  well-formed  longshore  bars  mirk- 
both  areas,  but  in  the  north  those  bars  are  stial lower  and  closer  to  shore  (ail 
within  tiie  first  400  meters).  South  ot  the  point  the  barred  zone  is  about  win 


Figure  7.  He.irh  an.l  nearshore  ti.it  hyn..  t  rv  In  t  lit*  vicinity  ot  the  I’entwater 
Channel  which  transects  the  heich  anil  barred  areas.  Hummocky 
bat  tivme  1 1  v  offshore  ir.iv  represent  material  dredged  from  the 
channel  ao.i  dumped  ottsliire  in  previous  vears. 

meters  wide  and,  as  in  the  north,  the  first  persistent  bar  lias  a  depth  ot  1  to 
-  meters  over  the  crest;  the  increase  in  depth  between  bars  in  an  offshore 
direction  is  taster  so  the  outermost,  well-developed  bar  has  5  to  b  meters 
over  the  crest  as  compared  to  i-  to  4-meter  depths  north  of  the  point.  The 
difference  in  shoreline  orientation  north  and  south  of  the  point  is  approxi¬ 
mately  Oil0,  which,  by  altering  the  nearshore  wave  conditions,  could  be  respon¬ 
sible  tor  the  contrast  in  bar  geometry  between  these  two  sections. 

a .  bar  Migration. 

> >n  the  basis  of  profile  changes  between  I9b7  and  1909  at  a  few  stations  In 
the  immediate  vicinity  of  the  He nt water  jetties,  Saylor  and  Hands  ( 1 9 7 1  • )  pro¬ 
posed  that  longshore  bars  migrate  landward  as  lake  levels  rise,  and  by  doing 
so  maintain  constant  depths  beneath  the  gradually  rising  lake  surface.  This 
proposal  was  consistent  with  Keulegan's  (194H)  conclusions  alter  studying  the 
tailors  controlling  har  t urinal  inn  in  wave  tanks.  However,  the  proposal  was 
contrary  to  a  1 1  the  other  tie  Id  studies  on  the  (.rent  Lakes.  Reports  ot 
previous  fieldwork  emphasized  long-term  stability  ot  bars  deeper  than  the 
ephemeral  bar  nearshore  (e.g.,  Davis  and  Mct.eary,  1901))  or  indicated  that 
during  periods  ot  persistently  rising  water  new  bars  were  continually  created 
Inshore  ot  the  old  series,  thereby  replacing  outer  bars  which  were  lett 
stranded  too  deep  to  be  allected  by  surface  wave  action  (tvans,  1940). 


Before  summarizing  the  bar  migration  observed  in  this  study,  consider 
that,  in  general,  any  interpretation  ot  profile  change  usually  entails  a 
belief  that  the  protile  occupied  only  those  positions  intermediate  between  the 
positions  determined  during  the  actual  survevs.  Whether  or  not  possible 


extreme  excursions  between  surveys  can  safely  be  ignored  depends  on  the  energy 
conditions  and  the  temporal  and  spatial  scales  involved,  as  well  as  on  the 

application  at  hand.  To  reveal  trends  in  long-term  bathymetric  adjustment  to 
higher  water  levels,  profiles  taken  even  several  years  apart  are  quite  useful, 
especially  if  all  the  profiLes  are  obtained  during  similar  phases  in  any 

seasonal  cycles  but  over  a  period  of  persistent  annual  change  (see  Fig.  8). 
The  magnitude  of  weekly  to  monthly  profile  changes  is  represented  in  Figure  8 
by  four  surveys  spanning  the  period  from  19  July  to  12  August  1967.  Their 
relatively  close  agreement  contrasts  with  the  difference  that  develops  between 
spring  and  fall  as  shown  by  the  May  and  August  1969  surveys.  The  long-term 

trend  in  bar  migration  can  be  seen  in  the  overall  change  from  19b7  to  1969  and 
in  the  comparison  of  bar  positions  in  those  years  with  the  final  bar  position 
determined  at  this  location  in  1975.  Careful  measurements  are  necessary  to 

discern  the  small  weekly  changes  from  possible  profile  error.  However,  the 
cumulative  effect  of  long-term  migration  clearly  exceeds  both  the  margin  of 
error  and  the  range  of  short-term  fluctuations. 


Figure  8.  Short-term  changes  In  bathymetry  versus  long-term  bar  migration.  The 
short-term  changes  are  Illustrated  by  four  surveys  In  1967.  May  and 
August  surveys  in  1969  reveal  larger  changes.  The  final  surveyeJ  bar 
position  In  1975  Illustrates  the  el  feet  of  long-term  bar  migration. 

Details  of  intermediate  surveys  at  other  stations  are  shown  in  Appendix 
A.  To  simplify  the  presentation  of  general  trends,  only  the  earliest  and 
latest  surveyed  bar  positions  are  shown  in  Figure  9.  The  original  survey  in 
1967  covered  the  area  in  the  immediate  vicinity  of  the  Pentwater  jetties 
(stations  3  to  8).  The  remaining  stations  (1,  2,  and  9  to  29),  spread  over 
the  adjacent  50  kilometers,  were  first  surveyed  in  19b9. 

Continued  monitoring  oi  profile  development  throughout  the  remainder  of 
the  rising  phase  in  lake  levels  and  for  several  years  thereafter  (until  1976) 
confirms  the  original  proposal  (Saylor  and  Hands,  1970)  that  bars  tend  to  rise 
with  the  water  level.  However,  landward  migration  of  the  bars  was  confirmed 
onLy  by  the  two  to  three  inner  bars  within  250  meters  of  shore.  The  outer 
bars  did  not  reveal  the  same  tendency  toward  shoreward  migration  as  the  inner 
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bars  did,  at  least  not  between  19b9  and  1975.  Furthermore,  the  outer  bar 
during  this  period  generally  lost  relief  as  a  result  of  various  combinations 
of  crest  erosion  (predominant  at  stations  2,  12,  16,  and  18)  and  trough  fill¬ 
ing  (predominant  at  stations  4,  8,  11,  and  19).  By  1975  these  two  processes 
had  progressed  to  the  point  of  completely  eliminating  the  outermost  bar  at  6 
of  the  33  stations  (4,  9,  11,  13,  19,  and  24). 

5.  Depths  of  Profile  Closure. 

All  profiles  collected  were  examined  for  evidence  of  a  limiting  depth 
below  which  there  were  no  bottom  changes  over  the  period  of  study.  In  19b7 
the  bathymetric  surveys  terminated  at  the  5-meter  contour.  Over  the  2-year 
period  from  1967  to  1969,  substantial  bottom  changes  occurred  throughout  the 
zone  from  the  shoreline  to  the  5-meter  contour.  This  evidence  of  deep  profile 
fluctuation  prompted  the  extension  of  surveys  to  greater  depths — 11  meters  in 
1969  and  21  meters  in  1971  and  1975.  Although  probable  depth  error  increased 
with  distance  from  shore,  the  longer  profiles  converged  at  their  outer  ends. 

Because  relief  on  the  longshore  bars  increased  significantly  from  one  bar 
to  the  next  in  the  lakeward  direction  and  the  bars  migrated  yearly,  the  enve¬ 
lope  of  bottom  change  also  increased  from  the  shoreline  lakeward  across  the 
barred  zone  (Fig.  10).  Beyond  the  outer  bar,  the  envelope  of  bottom  change 
narrowed  abruptly. 


Figure  10.  Envelope  of  bottom  change,  station  9  (  1967-75).  After  several  years  of  pro¬ 
file  adjustment  to  higher  water  levels,  the  envelope  of  changes  in  bottom 
elevation  is  thickest  in  the  zone  traversed  by  the  largest  migrating  bar  and 
narrows  abruptly  above  the  wave  uprush  and  below  the  barred  zone. 
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An  examination  of  all  profiles  indicated  that  instead  of  choosing  .1  si-igl,- 
limiting  depth,  it  would  be  more  realistic  to  pick,  two  depths:  the  critical 
depth,  a  shallower  depth  above  which  bottom  changes  typically  exceeded  <>.  i 
meter,  and  the  pinch-out  depth,  a  deeper  depth  below  which  there  is  no  evi¬ 
dence  of  change.  Between  the  critical  and  the  pinch-out  depths  there  is  a 
small  but  consistent  evidence  of  aggradation  (about  H.2  to  ''.1  meter  in  • 
years),  indicating  transport  and  accumulation  of  sediment  beyond  the  barred 
zone.  Beyond  the  pinch-out  depth,  changes  were  haphazard  and  genera  llv  less 
than  0.10  meter.  Which  of  the  two  indicators  of  closure  (critical  >r  pinch- 
out)  will  be  the  most  relevant  depends  on  the  application.  For  example,  when 
planning  a  sediment  budget  the  bathymetric  surveys  should  run  to  at  least  the 
pinch-out  depth.  Or,  the  other  hand,  when  selecting  a  site  for  placement  ol  a 
bottom  structure  or  instrument  package,  going  beyond  the  critical  depth  nay  he 
enough  to  preclude  burial  by  the  normal  processes  of  sedimentation. 

The  selection  of  the  closure  depths  involves  an  acknowledged  subjective 
evaluation  as  to  where  the  profiles  appear  to  close  at  each  survey  station. 
The  degree  to  which  individual  judgment  affects  these  estimates  is  illustrated 
in  Figure  11  by  two  estimates  obtained  independently  by  two  different  individ¬ 
uals  at  each  of  the  profile  stations.  Discrepancies  between  individuals, 
while  substantial  at  some  stations,  do  not  have  an  unacceptable  effect  on  the 
average  depth  for  a  broad  reach  of  shore.  Thus,  attempts  to  obtain  greater 
apparent  objectivity  in  the  selection  of  the  individual  depths  seem 
unnecessary . 


Stofion  Nos. 


Figure  11.  Depths  of  closure.  Estimates  of  profile  closure  were  made  inde¬ 
pendently  by  two  different  individuals  (represented  as  0  and  0) 
at  profile  stations  I  to  29.  Solid  lines  connect  their  estimates 
of  critical  and  plnch-out  depths  at  the  same  stations.  At  sta¬ 
tions  where  the  profile  did  not  extend  deep  enough  to  permit  a 
confident  selection  of  the  plnch-out  depth,  a  dotted  line  extends 
2  meters  below  che  critical  depth  estimate. 


21 


The  critical  depth  averaged  a  little  more  than  7  meters  on  the  north  side 
of  Little  Sable  Point,  and  a  little  more  than  8  meters  on  the  south  side 
(stations  16  to  29,  see  Fig.  9).  North  of  the  point,  the  average  pinch-out 
depth  was  10  meters.  South  of  the  point,  there  was  no  clear  pinch-out  on 
several  of  the  stations  because  the  1969  profiles  were  too  short;  where  the 
pinch-out  was  identified  it  averaged  11.5  meters.  Thus,  both  definitions 
suggest  deeper  profile  closure  south  of  the  point. 

The  individual  profiles  in  Appendix  A  may  be  useful  to  the  engineer  in 
determining  depths  of  measured  changes. 

6 .  Volume  Changes. 


a.  Stations .  To  test  the  assumption  that  the  volume  of  sand  eroded  from 
the  upper  beach  during  recession  was  matched  by  an  equal  volume  deposited  off¬ 
shore,  the  cross-sectional  area  between  profiles  at  each  station  was  cal¬ 
culated.  The  earlier  profiles  were  usually  too  short  to  include  all  of  the 
active  zone,  so  most  area  determinations  are  based  on  changes  between  1971  and 
1975.  Because  only  16  stations  were  reprofiled  in  1971,  4  of  the  longest  19b9 
profiles  were  used  to  supplement  the  area  change  measurements. 

b.  Calculations.  The  profiles  selected  for  volume  calculations  were 
digitized  at  5-meter  intervals  in  the  horizontal  from  the  landwardmost  to  the 
lakewardmost  points  common  to  both  the  earlier  and  later  surveys  (e.g.,  Fig. 
12).  The  results  of  all  the  volume  calculations  are  given  in  Appendix  B.  The 
difference  between  the  sequences  resulting  from  digitization  provides  a  se¬ 
quence  of  changes,  with  positive  values  indicating  a  fill  and  negative  values 
indicating  a  cut.  Multiplying  the  digitizing  interval  (5  meters)  times  the 
summation  of  all  elements  in  the  change  sequence  gives  a  measure  of  the  net 
volume  change  per  unit  width  alongshore.  If  the  elements  in  the  difference 
sequence  are  summed  from  their  landwardmost  point  to  some  arbitrary  point 
offshore,  the  product  of  that  sum  and  the  digitizing  interval  gives  the  net 
change  in  volume  per  unit  width  over  that  arbitrary  span.  Below  each  set  of 
digitized  elevations  there  is  a  continuous  curve  showing  the  change  in  volume 
per  unit  width  from  the  innermost  point  to  each  succeeding  point  across  the 
entire  active  profile.  This  cumulative  volumetric  curve  is  drawn  to  the  same 
horizontal  scale  as  the  profile. 

A  dashed  curve  plotted  on  the  same  axis  shows  the  average  thickness  of  the 
net  volume  change  if  it  were  distributed  uniformly  from  the  innermost  point  to 
the  end  point  for  which  the  change  was  summed  (Fig.  12). 

c.  Results .  Inevitably,  on  a  receding  shore  the  cumulative  volume  curve 
is  negative  from  the  inner  point  out  beyond  the  shoreline,  indicating  net 
degradation  or  cutting  over  the  upper  beach.  Small  zones  of  aggradation  or 
fill  offshore  cause  the  cumulative  volume  curve  to  increase  toward  zero  (Fig. 
12).  At  a  point  farther  offshore  the  cumulative  volume  curve  returns  to 
zero.  Between  this  balance  point  and  the  backshore  the  cut  and  fill  exactly 
balance  each  other;  i.e.,  neglecting  compaction  and  expansion,  the  sediment 
could  have  been  redistributed  within  that  zone  without  requiring  any  gain  or 
loss  to  the  outside.  Offshore  from  this  balance  point  the  cumulative  volume 
curve  would  ideally  not  depart  significantly  from  zero.  With  the  real  pro¬ 
files,  however,  the  cumulative  volume  curve  offshore  often  increases  about  as 
far  above  zero  as  it  was  below  zero  inside  the  first  balance  point  (App.  B). 
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Figure  12.  Schematic  diagram  of  volume  change  calculations. 

Farther  offshore  the  cumulative  volume  curve  usually  crosses  the  zero  line 
several  times  before  finally  smoothing  out.  A  positive  cumulative  volume 
curve  indicates  additional  sediment  was  supplied  from  outside  the  profile 
area.  If  this  sediment  came  from  offshore,  the  cumulative  volume  curve  would 
approach  zero  again  if  extended  to  the  pinch-out  depth.  As  it  turned  out, 
changes  in  depth  over  the  4-  to  b-year  period  were  so  small  near  the  pinch-out 
depth  that  the  total  volume  change  summed  over  the  entire  active  beach  was 
relatively  unaffected  by  deliberate  extensions  of  the  cumulative  volume  curve, 
and  therefore,  even  more  Insensitive  to  actual  uncertainties  encountered  in 
selecting  the  pinch-out  depth  (Fig.  11).  Although  the  cumulative  volume  curve 
approached  a  constant  near  the  pinch-out  depth,  it  usually  was  not  zero.  The 
value  of  the  cumulative  volume  curve  at  the  pinch-out  depth,  representing  the 
net  change  summed  over  the  active  profile,  is  tabulated  by  station  from  north 
to  south  in  Table  l.  The  concept  that  equivalent  volumes  are  eroded  from  the 
upper  beach  and  deposited  offshore  (sometimes  called  Bruun's  rule)  is  clearly 
invalid  when  applied  to  single  profiles.  In  fact,  witli  the  given  profile 
spacing,  there  is  no  sediment  balance  even  when  volume  changes  are  calculated 
over  reaches  several  kilometers  in  length.  Rosen  (1978)  pointed  out  similar 
local  imbalances  in  the  Chesapeake  Bay. 
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Table  1.  Volume  change  per  unit,  width  of  allure. 


Stat ion 

No . 

Time 

between 

surveys 

(vr) 

Volume  changes  (m^/m/yr) 

Zone  1 ^  Zone  2^  Zone  3^ 

1 

4.2 

-11.1 

-18.8 

-8.7 

2 

4.3 

-6.6 

-10.6 

2.3 

3 

4.3 

-11.7 

-15.5 

-20.4 

4 

4.3 

-6.8 

-9.4 

37.2 

5 

4.3 

2.6 

1.4 

37.6 

b 

4.3 

-0.5 

-1.9 

30.0 

7 

4.3 

-6.0 

-14.8 

-3.5 

8 

b.  2 

-4.8 

-9.0 

5.6 

10 

4.3 

-12.2 

-16.2 

9.6 

1  1 

5.8 

0.0 

-2.4 

26.9 

12 

6.3 

0.3 

-0.2 

12.0 

13 

4.2 

-16.5 

-17.9 

11.6 

18 

4.2 

-3.3 

-17.4 

18.1 

19 

4.2 

-11.5 

-20.3 

4.5 

20 

6.2 

-7.3 

-8.4 

3.7 

24 

4.2 

-0.9 

-8.0 

-2.4 

2b 

4.2 

-8.7 

-10.4 

10.4 

27 

5.9 

-10.2 

-13.1 

-44.3 

29 

4.2 

0.0 

-1.9 

-65.2 

Avg. 

-6.1 

-10.3 

3.4 

^Extends  from  dune  to  shore. 


Extends  from  shore  to  bottom  of  first  longshore 
trough. 

3 

Extends  over  entire  profile. 


For  the  entire  study  area,  the  total  volume  lost  from  the  upper1  beach  in  A 
years  (column  4,  Table  1)  averaged  41  cubic  meters  per  meter.  About  59  per¬ 
cent  of  this  loss  reflects  erosion  above  the  water  surface  (column  3,  Table 
1),  the  remaining  41  percent  occurs  between  the  shore  and  the  first  longshore 
bar.  The  net  volume  change  summed  over  the  entire  active  profile  produced  a 
small  positive,  but  statistically  insignificant  net  gain  of  sediment  (3  cubic 
meters  per  meter  of  shore  per  year).  Thus,  within  the  overall  survey  area, 
there  was  a  volumetric  balance  between  erosion  on  the  inner  part  of  the  beach 
and  deposition  offshore.  This  balance  suggests  that  future  profile  adjust¬ 
ments  to  different  changes  in  water  level  may  be  predicted  by  a  simple  geo¬ 
metric  model  discussed  in  the  next  section.  This  situation  will  not  exist  on 
many  open  ocean  coasts  where  eroded  beach  sands  are  transported  landward  by 
overwash  or  wind,  or  are  carried  into  inlets  and  deposited  on  ebb  or  flood 
tidal  shoals. 

IV.  PREDICTION  MODEL 

This  section  presents  an  idealized  concept  of  profile  adjustment,  dis¬ 
cusses  objections  and  difficulties  with  applying  the  concept,  and  shows  how 
these  difficulties  are  minimized  for  the  present  data  set.  Application  of  the 
concept  using  actual  measurements  is  followed  by  generalizations  and  specific 
guidance  on  applying  the  concept  to  other  areas  of  the  Great  Lakes. 
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As  described  by  Bruun  (  1962)  a  rise  in  the  mean  elevation  of  the  water 
surface  tends  to  shift  the  equilibrium  sand  profile  landward.  As  water  levels 
rise  erosion  prevails  on  the  upper  beach  and  the  shoreline  retreats.  Con¬ 
ceptually,  the  erosion  supplies  material  to  build  the  outer  part  of  the 
responding  profile  upward.  It  is  assumed  that  the  initial  profile  shape  is 
reestablished  farther  inland  and  at  a  distance  above  its  initial  position 
equai  to  the  change  in  water  level  z  as  depicted  in  Figure  13.  Thus,  the 
ultimate  retreat  of  the  profile  x  can  be  calculated  given  the  dimensions  of 
the  responding  profiles,  X  and  Z,  and  a  measure  of  the  stability  of  the 
shore-eroded  material  in  the  outer  zone,  R^. 


x 


«(MS8  <2> 

Z 


(1) 


where  sg  (z)  =1  if  z  >  0  (i.e.,  water  rising),  or  sg  (z)  =1  if  z  <  U  (i.e., 
water  faiiing. 


Figure  13.  Schematized  view  of  profile  adjustment 
as  two  rigid  translations. 


One  method  of  estimating  the  proportion  of  shore-eroded  materiai  that  will 
be  lost  is  to  use  the  textural  characteristics  of  the  active  beach  as  a  guide. 
If  the  newly  eroded  deposits  have  a  size  distribution  identical  to  that  of  the 
sediment  in  the  active  zone  prior  to  the  water  level  change,  then  only  insig¬ 
nificant  amounts  will  be  lost  through  selective  transport  processes  and 
Ra  —  l.  If  a  part  of  the  eroded  material  is  finer  than  the  overlying  native 

beach  it  may  be  carried  far  offshore  in  suspension  where  it  does  not  contrib¬ 
ute  to  the  building  of  a  new  profile.  In  which  case  RA  >  1  and  additional 
shore  erosion  must  compensate  for  the  loss.  Thus,  the  situation  is  simiiar  to 
the  problem  of  calculating  the  overfill  ratio  for  a  beach  nourishment  pro¬ 
ject.  Hobson  (1977)  explains  how  to  compute  to  evaluate  the  suitability 

of  borrow  material.  The  same  procedures  apply  here  except  the  "borrow  mate¬ 
rial"  characteristics  must  be  based  on  a  composite  sample  of  the  eroding 
section  of  the  shore,  i.e.,  the  upper  beach  in  the  case  of  increases  in  lake 
level  since  it  is  supplying  sediment  to  the  lower  part  of  the  adjusting  pro¬ 
file.  If  the  water  level  declines  the  lower  part  of  the  responding  profile  is 
eroded  to  supply  material  to  prograde  the  upper  profile.  In  this  case  the 

"native  material"  characteristics  must  be  based  on  a  composite  sample  of  the 

lower  profile  (i.e.,  the  zone  of  offshore  erosion).  In  either  case  the  native 
material  characteristics  must  be  based  on  a  composite  sample  of  the  entire 
responding  profile  from  the  limit  of  wave  uprush  to  the  point  of  profile 
closure. 

If  the  engineer  concludes,  without  specific  textural  data,  that  all  of  the 
shore-eroded  material  will  remain  in  the  zone  of  profile  adjustment,  then 

Ra  =  1.  If  the  engineer  estimates  by  other  methods  that  only  P  percent  of 
the  eroded  sediment  will  remain  in  the  active  zone  then  RA  =  100/P. 

Equation  (1)  with  R^  =  1  was  applied  to  sea  level  rise  on  Plorida  beaches 
by  Bruun  (1962),  and  in  this  context  is  often  referred  to  now  as  bruun's  rule. 
It  is  not  so  much  a  rule  in  any  formed  sense,  as  it  is  a  statement  of  a  fairly 
simple  concept  based  on  assumptions  which  had  been  used  by  many  early  coastal 
geomorphologists.  However,  explicit  applications  of  the  concept  prior  to 

Bruun  (1962)  are  unknown.  Although  references  to  the  concept  are  frequent,  it 
is  still  rarely  used  for  predictive  purposes. 


Civen  the  long-term  effect  of  rising  sea  levels  throughout  most  of  the 
Northern  Hemisphere,  it  may  be  wondered  why  the  sediment  balance  approach 
(Bruun's  rule)  has  not  been  more  widely  applied.  The  following  difficulties 
have  been  encountered  with  this  approach: 

(a)  Skepticism  as  to  the  adequacy  of  an  equilibrium  model  for 
explaining  short-term  dynamic  changes; 

(b)  difficulty  in  determining  RA  or  the  percentage  of  sediment 
lost  from  the  active  zone; 

(c)  problems  of  establishing  a  realistic  closure  depth  below 
which  water  level  changes  have  no  effect  on  profile  stability; 


(d)  confusion  arising  from  a  typographic  error  in  one  of  the 
equations  defining  profile  retreat  (Bruun,  1962);  and 

(e)  the  perplexity  caused  by  a  discontinuity  in  the  profile  at 
the  closure  depth  which  appeared  in  the  original  and  all  subsequent 
diagramatic  sketches  illustrating  the  concept. 

The  first  three  difficulties  (a,  b,  and  c)  warrant  serious  consideration 
before  applying  equation  (1);  items  (d)  and  (e),  although  perhaps  confusing, 
should  in  no  way  discourage  or  limit  use  of  equation  (1).  The  following 
paragraphs  address  each  of  these  difficulties  in  reverse  order. 

a.  Discontinuity  in  the  Profile  (Item  e).  Previous  diagrams  illustrate 
the  adjustment  of  a  profile  to  higher  water  levels  by  literally  disconnecting 
the  responding  part  of  the  bottom  from  the  static  region  offshore.  The  appar¬ 
ent  profile  discontinuity,  at  the  juncture  between  the  static  and  responding 
regions,  has  some  didactic  value  in  diagrams  to  the  extent  it  emphasizes  the 
congruency  between  initial  and  final  profile  shapes  in  the  active  region. 
Unfortunately,  it  also  creates  the  impression  that  the  model  is  inadequate  for 
explaining  the  transition  between  the  active  and  static  parts  of  the  profile. 
The  discontinuity  is  not,  however,  an  inherent  part  of  the  concept  but  rather 
an  artifice  of  the  diagrams.  Rigidly  translating  a  profile  upward  and  shore¬ 
ward  does  not  necessarily  lead  to  a  discontinuity  nor  even  a  change  in  slope 
as  is  demonstrated  later  in  this  report. 

b.  Error  in  an  Equation  (Item  d).  Bruun's  equation  (la)  (Bruun,  1962,  p. 
1 2 A )  is  dimensionally  incorrect  as  published.  This  error  may  have  discouraged 
some  readers  from  giving  Bruun’s  concept  their  full  consideration.  The  prob¬ 
lem  equation  is,  however,  unnecessary  to  the  development  of  this  concept 
(correctly  expressed  in  eq.  lb  of  Bruun,  1962).  The  validity  of  the  Bruun 
concept  and  of  equation  (1)  in  the  present  report  is  demonstrated  geometri¬ 
cally  in  Figure  13. 

Figure  13(a)  depicts  a  nearshore  profile  in  quasi-equilibrium  with  wave 
and  wave-related  forces.  Note  the  closure  depth  below  which  the  bottom 
presumably  does  not  adjust  to  surface  wave  and  current  conditions.  To  esti¬ 
mate  the  ultimate  shore  retreat,  the  adjustment  of  the  active  profile  is  then 
depicted  as  two  rigid  profile  translations. 

The  first  translation  moves  the  active  profile  (i.e.,  the  profile  between 
the  closure  depth  and  the  point  of  highest  wave  attack)  up  an  amount,  z,  and 
reestablishes  the  equilibrium  depths  below  the  elevated  water  surface  (Fig. 
13, b).  This  step  requires  a  volume  of  sediment  proportional  to  the  product 
of  X  (the  width  of  the  active  zone)  times  z  (change  in  water  level);  the 
volume  is  made  available  by  the  second  translation  which  is  recession  of  the 
profile  (Fig.  13, c).  Figure  13(c)  shows  that  x  units  of  recession  provide  a 
volume  of  sediment  proportional  to  the  product  of  x  times  2  (the  vertical 
extent  of  the  active  profile  from  the  critical  depth  up  to  the  average  eleva¬ 
tion  of  the  highest  erosion  on  the  backshore).  Equating  the  volumes  produced 
and  required  per  unit  length  of  shoreline  by  these  two  translations  (eq.  2) 
produces  equation  (1). 


then 


(Z) 


xZ  =  zX 


zX 


In  reality,  both  translations  occur  simultaneously  with  the  result  that 
the  closure  point  actually  migrates  upslope  as  the  water  level  rises.  Shift¬ 
ing  the  closure  point  upward  and  shoreward  will  affect  the  outcome  of  volumet¬ 
ric  calculations,  and  there  are  at  least  two  ways  to  account  for  this  small 
defect  in  the  geometric  justification  just  given  for  equation  (1). 

First,  a  closure  depth  midway  between  the  original  and  final  depths  could 
be  used  with  equation  (1)  to  improve  the  accuracy  of  the  calculation.  The 

horizontal  translation  of  the  profile  would  then  imply  a  slight  irregularity 

or  "step"  where  the  new  and  old  profile  shapes  meet.  The  step  would  consist 
of  a  wedge  of  surplus  sediment  above  the  closure  uepth  and  an  equal  volume 

deficiency  below;  therefore,  a  local  exchange  of  sediment  is  easily  imagined 
which  would  eliminate  the  step  and  completely  reestablish  the  identical  smooth 
profile  shape  without  affecting  the  overall  sediment  balance  expressed  in 

equation  (1).  This  method  of  accounting  for  the  migration  of  the  closure 
depth  is  easy  to  visualize  and  consistent  with  the  geometric  derivation  given 
for  the  predictive  equation. 

Second,  a  more  formal  development  of  the  sediment  balance  would  have 
integrated  between  profiles,  allowing  the  closure  point  to  move  in  infinites¬ 
imal  teps  with  the  water  surface.  This  approach  also  eliminates  the  step 
problem  and  results  in  the  more  precise  relationship: 

x  =  X  in  —  (3) 

L  -  Z 

Neither  method  of  adjusting  equation  (1)  (by  measuring  the  critical  depth 
from  an  intermediate  water  level  elevation  or  using  eq.  3)  is  generally  neces¬ 
sary  because  the  change  in  water  level,  z,  is  usually  so  small  relative  to 
the  total  height,  Z,  that  all  three  methods  provide  essentially  the  same 
results.  For  example,  if  z  <  G.1Z  all  results  agree  within  less  than  1 
percent . 

Thus,  the  simple  expression,  x  ~  zX/Z,  is  not  only  valuable  as  a  close 
approximation,  but  also  most  useful  because  it  is  easily  (a)  recalled  by 
visualizing  the  adjustment  of  two  rigid  translations,  (b)  expiained  in  the 
same  manner,  and  (c)  used  as  a  quick,  mental  check  on  the  ultimate  retreat 
expected  for  various  values  of  the  independent  variables. 

c.  A  Realistic  Closure  Depth  (Item  c).  Determining  a  realistic  closure 
depth  is  usually  extremely  difficult.  The  most  direct  ap,  ach  is  to  compare 
historic  bathymetric  surveys  of  the  site  in  question.  Unfortunately,  adequate 
survey  data  of  this  type  are  rare.  Neither  pier  nor  stadia  surveys  extern} 
deep  enough,  and  if  a  hydrographic  survey  does  extend  to  deep  water,  allow¬ 
ances  must  be  made  for  the  fact  that  both  sounding  errors  and  boat-positioning 
errors  usually  increase  significantly  with  depth  and  with  distance  from  shore. 
It  is  thus  often  impossible  to  substantiate  apparent  offshore  changes.  On  tin 


oceans,  waves  and  tides  create  di  t  t  ien  1 1  ies  in  cst.tbl  isliiap  a  datum.  \t  t  i  m  , 
it  is  impossible  to  distinguish  tin-  el  feet  or  1  .mg -pc  r  i  ud  swell  m  .  r>  ,n  n. 
waves  on  the  sounding  record  iron)  actual  sand  waves  .hi  the  sea  l  1  . . > r  C.'ui,  ...u 

and  Sarlin,  1  970;  Bruno  and  (lable,  1  ^  - 

Kinnllv,  it  deep  reliable  protiles  are  tumid,  evolution  it  tiie  1 . . 1 1 1  . 
depth  tor  use  in  equat  ion  (1)  requires  two  surveys  separated  by  an  ap,.1'1’!  :  1  ,4 
time  period  during  which  profile  adjustment  actual Iv  occurred  ia  response  to  , 
known  change  in  water  levels.  Hallerraeier  (  1  7  7 )  demuust  rates  t  ue  di;.e;,h  ..  . 

of  profile  closure  on  local  wave  conditions.  The  ditterence  helween  t  ,s  a 
closure  .it  two  sites  with  identical  wave  and  sediment  c  h  i  r.i.  t  e  r  i  s  t  l .  s  ,  .  •  u  with 
,1  stable  mean  water  level  and  the  other  with  a  recently  displaced  water  level , 
has  not  been  studied.  It  seems  plausible  that  storm  waves  .  unld  cans*  <  let 

profile  change  where  equilibrium  had  been  perturbed  by  toe  reient  shitt  in  t ' . 

mean  water  level,  and  yet  cause  onlv  sediment  notion  and  (.almost  by  dot  iuitiou) 
no  net  change  in  bottom  elevation  where  the  protile  was  in  equilibrium  with  a 
constant  water  level.  11  this  is  the  case,  real  water  level  changes  un¬ 
essential  it  repetitive  profiles  are  to  reveal  a  closure  depth  suitable  tor 

testing  the  Broun  concept.  (ilearlv,  many  problems  plague  the  delerminut  i-.n  ..1 
the  appropriate  closure  depth  and  theretore  discourage  .qq  1  i . at  i  on  ot  nr  aim's 
concept  for  predicting  future  shore  retreat. 

d.  Ka  or  the  Percentage  of  Sediment  Loss  (Item  bj .  squat  ion  (1)  .  m  be 

adjusted  to  account  for  any  sediment  lost  from  the  active  protile,  hut  onlv  i! 

the  volume  losses  can  be  determined.  Often  they  cannot.  Loss  occurs  when 
there  is  an  uncompensated  exchange  ot  sediment  beyond  the  surveyed  bound  tries. 
Losses  can  occur  offshore,  onshore,  or  alongshore.  i*n  the  w>est  -  oast,  subma¬ 
rine  canyons  complicate  the  determi n at  ion  ot  olt shore  losses.  On  the  ,yi  1  t 
coast,  hurricane  processes  have  moved  coarse  sediment  I  rota  as  deep  as  dm  meters 
onto  barrier  islands  (Hayes,  1 9  fc>  7 ) .  Return  currents  at  ter  hurricane  passage 

reportedly  spread  a  1-  to  J-cent imeter  layer  ot  beach  sand  over  homogenous  suds 
S  kilometers  t  rum  shore;  even  thicker  layers  of  nearshore  silts  ad  muds  re¬ 
portedly  moved  much  farther  gulf  ward  as  turbidity  currents  (Hayes,  1 9 7 ) . 

Onshore  losses  are  a  problem  on  the  east  coast.  High  tides  and  severe  storms 

transport  beach  sand  to  the  bay  side  ot  birriers  at  rates  ranging  t  rotti  more 

than  40  cubic  meters  per  meter  within  individual  over wash  deposits  during 
single  storms  to  about  1  cubic  meter  per  meter  lor  long  stretches  ol  shore 
yearly  (Schwartz,  1975).  The  engineer  must  consider  the  contribution  ot  these 
or  other  processes  to  sediment  losses  over  the  period  ot  his  study.  It  '/!'  is 
found  to  be  the  not  exchange  of  sediment  in  time,  T,  across  the  boundaries  ot 
a  control  area  with  longshore  length,  Y,  then  the  anticipated  retre.it  should 
be  reduced  by  t)T/YZ: 

u>  .it 

x  =  - - - -  (4  ) 

Z  YZ 


e.  Adequacy  of  an  Equilibrium  Model  (Item  a).  l-se  of  equilibrium  assump¬ 
tions  to  model  dynamic  coastal  changes  also  deserves  scrutiny.  The  idea  ot  an 
"equilibrium  beach  profile"  has  had  a  long  history  (e.g.,  Kenneman,  1902);  how¬ 
ever,  opinions  still  differ  as  to  exactly  what  the  concept  actually  entails. 
By  one  definition,  the  profile  of  equilibrium  is  the  ultimate  shape  which 
coastal  processes  strive  to  impart  to  a  beach.  Of  coarse,  nature  seldom  re¬ 
mains  constant  long  enough  for  a  strict  equilibrium  to  develop.  In  the  present 
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context,  the  term,  r>(/n  r  '  mt  rlr‘<  !  e  refers  to  ,1  curve  oi  tixed  si/.e  aiui 
shape  which  "adequately"  represents  the  "average"  protile  shape  belore  per¬ 
turbation  by  a  shitt  in  water  level.  By  assumption,  shore  erosion  eventually 
returns  the  profile  to  this  same  shape  alter  it  is  displaced  as  a  result  ut 
the  water  level  change  (see  Fig.  11). 

v  willingness  to  accept  equilibrium  as  a  reasonable  approximation  is  not 
inconsistent  with  recognition  of  seisonal,  storm,  or  other  temporary  profile 
fluctuations.  Careful  judgment  should  be  made  on  a  case-bv-case  basis,  il 
tield  protiles  claim  to  represent  quas  i  -e<ju  i  1  i  br  i  uin  conditions.  l.eueralls, 
the  claim  will  he  more  reasonable  the  longer  the  time  t raine  ot  the  study.  The 
spatial  extent  of  the  study  is  also  important.  Usually,  the  longer  the 
stretch  ot  shore,  the  more  likely  that  longshore  variations  will  also  average 
out,  thereby,  providing  an  overall  equilibrium. 

This  discussion  has  shown  that  the  Hruun  concept  is  theoret  ica  l  Iv  sound 
but  difficult  to  apply  in  ttie  tield.  The  next  subsection  examines  how  some  ot 
the  difficulties  discussed  above  are  avoided  in  the  lake  Michigan  data. 

1 .  Suitability  of  Present  Data  for  Testing  the  Sediment  Balance  Concept. 

The  ways  that  previously  discussed  difficulties  (items  a  to  e)  ultect  Like 
Michigan  data  are  outlined  here,  before  an  actual  application  ot  the  data  in 
tlie  next  subsection.  Difficulties  (d)  and  (e)  should  not  limit  ipplication  ot 
the  model  to  any  data  set  tor  reasons  discussed  in  the  Iasi  subsection. 

Establishing  a  realistic  closure  depth  (item  c)  depends  on  uccnrite 
repetitive  profiling.  Profile  errors  increase  with  distance  from  shore. 
Fortunately,  the  bottom  drops  oft  to  suitable  depths  relatively  rapidly  in  the 
present  study  area.  Furthermore,  the  llre.it  Lakes  are  tree  t  rom  l  i  da  1  varia¬ 
tions  as  we  1 L  as  from  long-period  swell.  The  Lreat  Lakes  are  notorious  for 
their  large  storm  surges  and  seiches;  however,  based  on  extensive  water  level 
measurements  in  1969,  it  was  concluded  that  these  disturbances  are  not  a 
signifi  ant  problem  in  the  present  study.  By  choosing  the  right  time  ot  year 
and  surveying  only  when  conditions  are  calm,  it  Is  possible  to  avoid  datum  ind 
bottom  ambiguities.  Note  the  absence  of  confusing  wave  interference  on  tin- 
raw  fathogram  in  Figure  14. 

The  difficulty  of  determining  sediment  losses  (item  b)  on  the  Lre.it  Lakes 
is  greatly  simplified  by  the  absence  of  submarine  canyons,  hurricanes,  and 
overwash  events.  Fluvial  sediment  input  is  also  no  problem  because  all  rivers 
entering  eastern  Lake  Michigan  flow  through  deep  inland  sediment  traps. 
Dredging  at  Pentwater  ChanneL  is  well  documented.  On  the  average,  bO,OUU 
cubic  meters  is  removed  annually,  and  some  of  this  is  returned  to  adjacent 
beaches.  Inlet  losses  have  only  a  small  effect  on  the  overall  sediment  budget 
for  the  broad  study  area.  Thus,  in  the  present  application  Q  (eq.  4)  will 
have  a  negligible  effect. 

The  only  process  supplying  new  sediment  to  the  active  profile  is  shore 
recession.  Furthermore,  shore  deposits  and  backshore  bluffs  within  the  study 
area  contain  less  than  1  percent  silt,  making  it  unnecessary  to  correct  for 
any  unstable  fine  fraction  (i.e.,  =  1;  eq.  1).  Thus,  a  number  of  site- 

specific  attributes  simplify  sediment  balance  for  the  study  area. 
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Tne  dependence  ut  the  model  on  equilibrium  assumptions  (item  a)  makes  it 
difficult  to  get  convincing  field  confir’m.t-'rm  because  waves,  currents,  and 
conditions  of  sediment  supply  never  remain  constant  for  long.  on  tiie  other 
hand,  because  it  is  an  equilibrium  model  the  potential  tor  irq  '  ’  ; '  ’  ’c  is 

broadened.  The  model  can  even  apply  to  situations  where  storms  or  nun’s 
influence  upset  equilibrium.  The  predicted  retreat  in  such  cases  would 
indicate  the  adjustment  required  by  just  the  change  in  water  level;  eltects 
due  to  other  changes  would  have  to  be  super i mposed  it  signiticanl.  Thus,  once 
the  concept  is  confirmed,  equilibrium  tends  to  make  application  easier. 

The  model  itself  provides  no  indication  ot  the  time  period  required  lor 
the  beach  to  return  to  equilibrium.  Errors  in  misjudging  equilibrium,  and 
failing  to  account  for  the  lag  between  cause  and  effect,  are  all  too  easy  to 
make  if  the  data  cover  only  a  small  reach  of  shore  or  a  short  period  ot  time. 
The  length  of  time  and  the  number  of  profiles  studied  here  are  thought  tu  be 
sufficient  to  avoid  this  problem. 

4.  Application  of  the  Sediment  balance  Approach. 

a.  Longshore  Contributions.  Wave  data  suitable  tor  prediction  til  long¬ 
shore  transport  rates  are  not  available  in  the  study  area.  Various  indica¬ 
tions  of  the  direction  of  transport  are  compiled  in  Figure  11.  f.vidence  Iron, 
coastal  geomorphology  (Hands,  1970),  from  longshore  changes  in  grain  size 
(Saylor  and  Hands,  1970),  from  the  pattern  of  channel  shoaling  (Hands,  197ba), 
and  from  data  hindcast  for  extreme  storms  (Resio  and  Vincent,  197bc)  suggests 
that  the  direction  of  longshore  transport  in  the  vicinity  ot  Pentwater  Harbor 
is  predominantly  southward,  but  subject  to  frequent  reversals;  extrapolation 
from  SavilLe’s  (1953)  hindcast  data  suggests  a  northward  transport.  Littoral 
Environment  Observation  (LEO)  data  from  Mears  State  Park  were  inconclusive — 
too  short  a  record  and  subject  to  the  effects  of  a  large  eddy  and  reflected 
waves  from  the  Pentwater  jetties.  Near  profile  station  17,  the  extreme  storm 
data  and  the  usual  deflection  of  Silver  Lake  Creek  crossing  the  beach  suggest 
that  the  direction  of  transport  changes  to  northward  on  the  south  side  ot 
Little  Sable  Point.  Beyond  the  southern  limit  of  the  study  area,  storm  data 
from  White  Lake  and  Muskegon  suggest  a  close  balance  between  northward  and 
southward  flows  in  that  region.  South  of  brand  Haven  the  goomorphol ogy  and 
storm  data  indicate  net  southward  transport  for  the  remainder  oi  the  eastern 
shore.  Therefore,  there  is  a  consistent  pattern  ot  drift  moving  toward  Little 
Sable  Point  from  the  north  (Summit  Park)  and  from  the  south  (White  Lake)  (Fig. 
15). 

Long-term  convergence  of  drift  toward  the  Silver  Lake  dunes  would  be 
consistent  with  the  evolution  of  Little  Sable  Point  from  a  shallow  embayment 
several  thousand  years  ago  when  water  levels  were  7  meters  above  modern  levels 
(Hough,  1958)  to  the  dune-covered  coastal  promontory  of  today. 

The  areas  from  Ludington  to  Summit  Park  and  from  White  Like  to  Muskegon 
appear  to  be  natural  boundary  zones  of  longshore  divergence  (Fig.  15).  In 
addition  to  these  natural  boundaries,  the  jetties  at  Ludington  and  at  the 
pumped  storage  facility  4  kilometers  farther  south  (Fig.  lb)  are  also  obsta¬ 
cles  to  sediment  input  from  the  north.  The  jetties  and  entrance  channel  .it 
White  Lake  likewise  reinforce  the  natural  southern  boundary.  Present-day 
processes,  storm  patterns,  and  engineering  projects  thus  limit  tla-  possible 
sources  of  drift  converging  toward  the  Silver  Lake  dunes  to  those  beaches  and 
bluffs  primarily  within  the  present  study  area  (see  Fig.  4). 
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figure  1  'i .  Direction  of  net  longshore  transport  along  the  cast.  rn  shore.  Nut 
all  sources  are  in  agreement,  hut  most  support  the  hypothesis  t  li.at 
there  are  areas  of  longshore  divergence  hevond  the  studv  irea  Noth 
to  the  north  and  south,  and  that  longshore  transport  within  the 
study  area  converges  near  l.i  t  I  le  Sahle  Point. 


Summit  Park  Pentwater 


figure  lu.  Pumped  storage  facility  south  or  Ludington.  The  shore  soutti 
ot  Ludington  is  a  zone  or  general  longshore  divergence.  The 
jetties  built  at  tire  pumped  storage  tacility  in  1  SI 7 1  turther 
restrict  longshore  transport  in  this  area  or  divergence  and 
Llius  establish  a  specific  northern  limit  tor  sources  ot  sand 
to  the  present  study  area. 

The  assumption  that  there  is  no  significant  longshore  input  or  Losses 

beyond  the  present  study  limits  is  reasonable,  especially  considering  the 

minor  impact  any  imbalance  would  have  on  the  5U-kiiometer  stretch  ot  shore 
during  this  period  of  rapid  shore  erosion.  Por  example,  a  net  inflow  ol 

100,000  cubic  meters  per  year  (an  improbably  large  tigure)  would  be  volu- 
metrically  equivalent  to  recession  of  only  0.J  meter  per  year,  (10  u  /  y  t* )  / 
(i  x  10"*  m  x  lU  tn),  while  the  observed  shore  recession  actually  averages 
J .  b  meters  per  year.  So  the  irviximuin  conceivable  longshore  input  is  small 
relative  to  the  enormous  exchange  of  sediment  onshore  and  offshore  during  this 
period  of  rapid  profile  adjustment. 

b.  A  Possible  Inland  Loss.  A  possible  inland  loss  on  Little  Sable  Point 
complicates  the  otherwise  simple  sediment  balance  lor  this  area.  The  Silver 
Lake  dunes  occupy  about  b  kilometers  ot  shoreline  between  profile  stations  li 
and  17  and  extend  more  titan  a  kilometer  inland  (Pig.  9).  These  actively 
migrating  dunes  reach  heights  ot  V>  meters  along  t In'  inland  halt  ot  the  dune 
field  (Pig.  17).  Along  the  shoreline,  the  dune  ridges  >  rest  about  7  meters 
above  lake  level  (Pig.  1H)  and  some  ponded  inlerdune  areas  are  at  approxi¬ 

mately  the  same  level  as  Lake  Michigan  (177  meters.  International  l.reat  Lakes 
Datum).  Por  ninny  hundreds  of  years  this  dune  Held  has  been  ted  by  the  con¬ 
vergence  of  longshore  transport  toward  Little  Sable  Point  t rom  both  the  north 
and  south,  and  by  the  inland  transport  ot  sand  by  prevailing  wist  winds.  The 

dune  shoreline  receded  more  than  any  ot  the  adjoining  beaches  during  the  study 

period  (Pig.  19).  The  vastness  ot  the  dune  field,  the  etteet  ot  I-  to  J-meter 

wave-cut  bluffs  which  marked  much  ot  the  dune  shoreline,  and  the  virtual 
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figure  1/.  Silver  Lake  dunes,  looking  from  Lake  Michigan  across 
the  dunes  to  Silver  Lake.  Spit  extension  probably 
sealed  off  the  two  bodies  of  water  soon  alter  the 
Nipissing  high  lake  levels  which,  according  to  hough 
(l9Bb),  would  have  been  about  J.UUU  years  oetore 
Present  (B.P.).  Continued  longshore  transport  and 
prevailing  west  winds  built  the  dune  field  along  the 
front  of  rhis  low  receptive  embavment. 

absence  of  an  exposed  beach  across  which  the  wind  could  blow  make  it  difficult 
to  estimate  the  volume  of  sand  actually  transferred  inland  during  the  recent 
period  of  high  water.  It  is  assumed  that  inland  losses  to  the  Silver  Lake 
dunes  hetween  1969  and  1976  exceeded  the  gain  of  sediment  supplied  to  the 
adjusting  profile  from  the  dunes  by  only  a  small  amount  which  can  be  neg lee  *  j 
in  the  calculation  of  an  overall  sediment  balance  for  the  larger  study  area. 

c.  Measured  and  Predicted  Shore  Retreat.  Because  the  initial  1 9 f> 7  survey 
covered  only  a  small  area  in  the  vicinity  of  Pentwater  Harbor,  an  area  subject 
to  less  recession  than  the  surrounding  "undisturbed"  beaches  (Hands,  1979), 
testing  of  the  sediment  balance  approach  was  best  done  hv  using  the  1969-76 
survey  data  and  excluding  measurements  made  within  BOO  meters  of  the  Pentwater 
jetties.  The  extent  of  shore  covered  (21  profiles  spread  over  BO  kilometers) 
and  the  length  of  time  monitored  (7  years),  together  with  the  sizable  increase 
in  mean  water  level  during  the  study  period  and  the  generally  near-ideal  con¬ 
ditions  discussed  previously,  make  this  application  the  most  realistic  field 
test  of  equilibrium  profile  migration  to  date. 

Measurements  of  the  width  of  each  profile  from  the  vegetated  dune  line  to 
the  pinch-out  depth  for  each  station  were  taken  and  averaged  to  obtain  X  =  921 
meters.  The  heights  of  ttie  scarps  which  waves  had  cut  in  t  lie  toredune  were 
also  measured  at  each  station.  As  the  profiles  had  not  vet  developed  a 
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Figure  19.  Longshore  variation  in  net  shore  retreat  from 
1969  to  1976. 

significant  scarp  in  1969,  the  average  dune  height  for  the  period  of 
adjustment  was  approximately  one-half  the  average  scarp  height  in  1976. 
Average  dune  height  added  to  the  average  pinch-out  depth  taken  from  Figure  11 
established  the  vertical  dimension  of  the  adjusting  profile  (the  X  in  eq.  1 
and  Fig.  13).  The  resulting  average  value  of  Z  was  13. b  meters.  Under  the 
discussed  assumption  of  regional  sediment  balance,  the  ratio  X/Z  times  the 
measured  water  level  change  ( z  =  0.20  meter)  equals  the  ultimate  shore  retreat 
(-13.6  meters).  The  retreat  actually  measured  between  1969  and  197b  also 
averaged  -13.6  meters.  Considering  the  measurement  and  sampling  errors  in¬ 
volved  in  determining  each  independent  variable  a  predictive  capability  of 
less  than  a  tenth  of  a  meter  certainly  is  not  claimed,  but  the  results  clearly 
confirm  the  appropriateness  of  the  equilibrium-sediment  balance  approach  when 
applied  in  the  proper  setting. 

As  noted  previously,  pinch-out  depths  are  deeper  south  of  Little  Sable 
Point  than  to  the  north.  The  eroding  dunes  are  also  higher  there,  which  even 
further  enlarges  the  vertical  dimension  of  profile  adjustment  south  of  Little 
Sable  Point.  Consequently,  the  equilibrium  prediction  might  be  applied  sepa¬ 
rately  to  the  two  regions.  Likewise,  because  additional  surveys  wore  con¬ 
ducted  in  1971  and  1975,  separate  predictions  could  be  applied  to  these 
shorter  time  intervals  (1969  to  1971  and  1969  to  1973)  as  well.  Thus,  par¬ 
titioning  the  original  data  provides  nine  individual,  though  not  independent 
tests  (Table  2).  The  greater  pinch-out  depth  south  of  Little  Sable  Point 
increases  both  the  width,  X,  and  height,  Z,  estimates  in  a  compensating 
fashion  so  that  there  is  little  effect  on  the  predicted  outcomes.  The  values 
predicted  for  north  of  Little  Sable  Point  are  essentially  the  same  as  pre¬ 
dicted  for  south  of  Little  Sable  Point  for  eacli  of  the  three  time  periods 
(Fig.  20).  Considering  prediction  versus  measurement,  the  predicted  retreat 
from  1969  to  1971  was  too  high  for  all  three  areas  (117  percent  high  lor  tin- 
area  as  a  whole).  The  prediction  for  19b9  to  1975  was  also  high,  hut  not  as 
far  off  as  before  (45  percent  high  for  the  whole  area).  These  ove res l i mu L  i .  ms 
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Profile  Retreol  (m) 


Table  2.  Predicted  and  observed  profile  retreat. 


Study  area 

Survey  periods 

1969-1971  1969-1975  1969-1976 

x  “  0.12  m  x  ■  0.39  m  x  ■  0.20  m 

X  -  870  m  X  =  1,020  m  X  =  923  m 

Average  height,  Z  (m) 

Northern  section 
(stations  1  to  15) 

Southern  section 
(stations  16  to  29) 

Whole  area 
(stations  1  to  29) 

10.84  12.15  12.50 

12.90  14.28  14.80 

11.86  13.21  13.60 

Predicted  retreat,  Xz/Z  (m) 

Northern  secton 
Southern  section 

Whole  area 

9.63  27.93  13.92 

9.49  27.86  13.78 

9.34  27.25  13.57 

Observed  retreat  x  (m) 

Northern  section 
Southern  section 

Whole  area 

4.6  20.0  12. b 

3.b  16.8  14.8 

4.3  18.8  13.6 

Overprediction  (pet) 

Northern  section 
Southern  section 

Whole  area 

109  40  10 

164  h6  7 

117  45  0 

30r 


( yr) 

Figure  20.  Calculated  versus  measured  retreat.  The  predicted  ultimate  retreat,  in 
response  to  post-1969  changes  in  mean  lake  level,  exceeded  the  observed 
retreat  by  more  than  100  percent  in  1971  and  about  50  percent  in  1975, 
presumably  because  the  active  profile  had  not  had  time  to  completely  read¬ 
just  to  the  higher  water  levels.  Almost  perfect  agreement  had  developed  by 
the  time  of  the  last  survey,  3  years  after  the  lake  levels  peaked. 
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of  retreat  are  attributed  to  the  fact  that  profile  retreat  was  actually  lag¬ 
ging  behind  the  lake  level  rise.  As  hypothesized  earlier  (Hands,  1976a), 
rising  water  levels  establish  a  potential  for  erosion  and  realization  of  that 
potential  requires  sediment  redistribution,  i.e.,  work  which  d-nends  on  the 
energy  being  available.  The  convergence  of  measured  and  predicted  retreat  in 
both  regions,  3  years  after  annual  lake  levels  had  stabilized,  suggests  that 
several  storm  seasons  may  be  required  to  readjust  the  profile  to  changes  in 
mean  water  level  of  several  tenths  of  a  meter. 


According  to  the  model,  which  works  well  here,  the  problem  of  predicting 
the  effect  of  lake  level  changes  is  equivalent  to  the  problem  of  identifying 
the  pinch-out  depth.  The  remarkable  confirmation  of  theory  and  data  in  the 
present  case  highlights  the  need  to  generalize  a  method  applicable  to  similar 
regional,  long-term  settings  but  where  wave  energies  and  therefore  pinch-out 
depths  might  be  significantly  different. 


5.  Using  Wave  Climate  to  Estimate  the  Pinch-Out  Depth. 

In  the  model,  the  closure  depth  is  the  point  below  which  the  bottom  does 
not  adjust  to  changes  in  water  surface  elevation.  In  the  field,  this  point 
was  approximated  by  averaging  the  upper  bounds  of  the  region  of  negligible 
profile  change  in  repetitive  surveys  (pinch-out  depths).  The  closure  depth, 
thus  established,  is  not  necessarily  appropriate  for  other  areas  of  the  Great 
Lakes.  The  depth  of  profile  closure  should  vary  regionally  with  the  wave 
climate.  Unfortunately,  the  repetitive  profile  record  is  usually  not  suffi¬ 
cient  to  establish  this  parameter. 

In  these  cases,  knowledge  of  the  wave  climate  is  useful.  Wave  gage  data 
obtained  during  profile  survey  periods  are  too  short  to  be  indicative  of  the 
important  long-term  conditions  in  the  study  area;  however,  wave  climate  data 
are  available  from  other  sources  including  hindcast  studies  (Saville,  1953; 
Resio  and  Vincent,  1976a,  1976b,  1976c,  1977,  1978),  shipboard  observations 

(Pore,  et  al.,  1971;  National  Oceanic  and  Atmospheric  Administration,  1975), 
U.S.  Coast  Guard  reports  (Liu  and  Housley,  1969),  and  the  LEO  program  (Weggel, 
1979).  Considering  site  specificity,  long-term  coverage,  and  the  availability 
of  comparable  data  for  the  entire  U.S.  shoreline  of  the  Great  lakes,  Resio  and 
Vincent's  reports  were  chosen  as  the  basic  reference  for  extrapolating  profile 
response  from  the  present  study  area  to  those  with  significantly  different 
wave  environments.  Their  wave  climate  parameters  were  generated  by  a 
numerical  hindcast  model  using  wind  data  from  the  extreme  storms  recorded  over 
a  30-year  period.  The  parameters  thus  describe  only  the  deepwater  storm 
conditions.  Because  the  maximum  depth  of  profile  response  depends  on  the 

higher  waves  and  because  only  a  consistent,  relative  measure  of  spatial  wave 

variability  is  needed,  the  milder  waves  though  important  in  profile 

development  need  not  be  considered  here.  It  is  reasonable  to  assume  that  the 
maximum  depth  of  intense  bottom  agitation  depends  on  at  least  the  wave  period 
and  the  shoaled  and  refracted  wave  height,  but  Hallermeier  (1977)  found  that 
the  maximum  depth  in  a  number  of  actual  design  wave  conditions  was  essentially 
proportional  to  deepwater  wave  height  alone. 

The  wave  height  data  necessary  to  estimate  the  pinch-out  depth  for  any 

Great  Lakes  site  are  given  in  Appendix  C.  The  average  pinch-out  depth 
established  within  the  present  study  area  is  2.1  times  the  average  5-year 


return-period  wave  height  for  the  area.  In  the  absence  of  direct  profile 
change  measurements,  the  pinch-out  depth  for  other  regions  is  thus  estimated 
HS  I  ^ 

d  ~  2.1  h5  (5) 

where  h^  is  the  5-year  return-period  height  given  in  Appendix  C. 

The  accuracy  of  this  approach  is  not  known.  For  an  idea  of  how  sensitive 
the  prediction  of  profile  response  is  to  errors  in  estimated  closure,  assume 
that  this  study  is  restricted  to  the  area  north  of  Little  Sable  Point  (sta¬ 
tions  1  to  15,  Fig.  9).  North  of  the  point,  the  pinch-out  depth  (from  direct 
measurements)  averages  10  meters.  With  no  repetitive  profiles  south  of  the 
point  the  estimate  would  be  that  pinch-out  occurs  at  1.96  h^.  The  1.96  would 
be  a  less  reliable  estimate  of  the  coefficient  in  equation  (5).  The  average 
h^  south  of  Little  Sable  Point  was  5.39  meters  (Michigan  stations  15  and  16, 
see  App.  C),  so  the  estimated  closure  depth  would  have  been  10.56  meters  which 
is  1.43  meters  or  12  percent  too  small  (Table  3).  Adding  this  value  to  the 
mean  dune  height,  D,  south  of  Little  Sable  Point  produces  a  new  estimate 
of  2;  the  distance  of  the  10.56-meter  contour  from  shore  produces  a  new 
estimate  of  X  (Table  3).  Using  these  new  values,  the  estimated  response  to 
a  0.2-meter  increase  in  lake  level  would  be  13.9  meters  which  is  only  1 
percent  over  the  value  obtained  from  actual  measurements  south  of  the  point. 
This  exercise  illustrates  the  self-compensating  tendency  which  errors  in  the 
pinch-out  depth  have  on  equation  (1).  The  depth  estimate  was  12  percent  too 
small  when  the  procedure  was  applied  to  data  different  from  those  used  to 
estimate  the  coefficient  k  in  d  =  kh^.  The  effect,  however,  was  to 
introduce  less  than  1  percent  error  in  the  predicted  shore  retreat. 


Table  3.  Cross-validation  Indicates  the  effect  of  estimating  plnch-out  depth  from  vave  climate 
_ data.  Calculated  estimates  (hatted)  are  compared  to  measured  values  (nonhatted). _ 


Z  -  d  +  D 


Whole  area  5.3  2.1  11 

(stations  1  to  29) 

Northern  section  5.1  1.96  10 

(stations  1  to  15) 


Southern  section  5.4  -  12  10.56  14.8  13.36  1,020  928  13.78  13.89 

(stations  16  to  29)^ 

Error - - - - - 12%  -  10%  -  91  -  <11 

*Used  for  an  Independent  determination  of  k. 

2 

Used  for  error  check  by  comparison  with  estimates  based  on  data  from  the  northern  section  only. 


To  be  realistically  applied,  the  model  should  have  input  from  many 
profiles  spaced  along  a  section  of  coast;  as  a  consequence,  there  is  little 
point  in  partitioning  the  present  data  set  any  further.  The  cross  validation 
shown  above  does  not  reflect  all  the  drawbacks  of  estimating  the  pinch-out 
depth  from  wave  climate  because,  using  adjacent  sections  of  coast,  it  does  not 
introduce  the  full  range  of  bathymetric  variability  nor  the  range  of  wave 
environments  within  the  lakes.  How  well  wave  climate  estimates  from  widely 
different  environments  will  perform  remains  uncertain.  However,  the  prospects 
seem  good  and  alternatives  nonexistent.  The  model  .  should  be  applied  cau¬ 
tiously,  and  wherever  there  is  any  indication  of  how  well  or  poorly  it  worked, 
the  results  should  be  reported. 
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6.  Inferences  from  ProfiLe  Shape  Alone 


a.  Preliminary  Check.  Prior  to  any  detailed  evaluation  of  the  variables 
required  as  input  to  the  model,  a  preliminary  examination  of  the  profile  shape 
near  the  presumed  closure  depth,  d,  will  indicate  how  reliable  these  evalua¬ 
tions  need  to  be.  If  the  profile  shape  changes  abruptly  near  this  depth,  then 
the  choice  of  d  may  strongly  affect  the  resulting  prediction,  depending  on 
whether  the  exact  value  chosen  is  above  or  below  the  break  in  slope  (see  Fig. 
21, a).  If,  however,  the  bottom  is  planar  and  sloping  at  the  right  angle,  the 
ratio  of  Z/X  and  therefore  the  predicted  retreat  will  be  unaffected  by  vari¬ 
ations  in  d  over  a  wide  range.  In  such  cases  the  exact  value  of  d  used  to 
evaluate  equation  (2)  will  be  unimportant  (Fig.  21, b). 


o.  UNMATCHEO  SLOPES 


b.  UNIFORM  SLOPE  a'  =  a> 


Figure  21.  Importance  of  offshore  slope.  In  case  a  the  prediction 
of  retreat  will  be  much  more  sensitive  to  the  correct 
selection  of  a  pinch-out  depth  than  in  case  b. 


b.  Is  the  Prediction  Conservative?  The  evaluation  of  d  on  the  basis  of 
wave  data  may  result  in  an  estimate  of  closure  which  the  engineer  feels  is 
either  too  low  or  too  high.  Yet  he  may  have  no  specific  evidence  on  which  to 
base  another  choice.  The  engineer  should  determine  if  the  suspected  error 
strengthens  or  weakens  arguments  based  on  the  sediment  balance  model.  Exam¬ 
ination  of  profile  shape  resolves  the  possibilities  as  shown  in  Figure  22. 
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Figure  22.  Diagram  for  determining  if  a  suspected  error  weakens  or 

strengthens  arguments  based  on  the  sediment  balance  prediction. 

If  a  line  connecting  two  comtemplated  closure  depths  extends  below  the  average 
height  of  the  dunes  throughout  the  section  of  shore  under  study,  then  over¬ 
estimating  the  closure  depth  causes  equation  (1)  to  overpredict  the  response; 
underestimating  the  closure  causes  equation  (1)  to  underpredict  the  response. 
But,  if  the  line  extends  above  the  dune  height,  then  overestimating  the 
closure  underpredicts  response  and  underestimating  the  closure  overpredicts 
response.  If  the  extended  line  intersects  the  dune  crest,  the  prediction 
remains  unaffected  by  the  error  in  d.  These  relationsips  will  apply 
regardless  of  whether  the  predicted  response  is  a  retreat  or  an  advance  of  the 
shore. 


c.  Ideal  Long-Term  Development.  The  longshore  variability  in  slope  near 
a  depth  of  d  is  another  item  to  check.  If  all  profiles  have  the  same  off¬ 
shore  slope,  a,  and  if  it  is  assumed  that  long-term  recession  unearthed 
deposits  similar  to  the  modern  substrate  and  the  wave  climate  has  been  sta¬ 
tionary,  then  tan  a  =  Z/X;  i.e.,  a  is  not  only  the  actual  bottom  slope,  but 
also  Hick's  (1972)  effective  angle  of  shore  response.  Profile  migration  in 
response  to  rising  water  levels  under  these  conditions  would  ideally  leave  a 
slope  below  its  trailing  edge  which  could  serve  as  a  clue  to  past  recession. 
If  the  shore  formerly  supplied  a  greater  volume  of  littoral  material  per  unit 
of  recession  (because  dunes  or  bluffs  were  higher  or  contained  a  larger  per¬ 
centage  of  suitable  littoral  material),  then  the  slope  beyond  the  trailing 
edge  would  tend  to  be  convex.  An  increasing  supply  of  sand  would  tend  to 
produce  a  concave  slope  as  erosion  provides  more  and  more  sediment  for  each 
unit  of  recession. 
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To  simplify  profile  representation,  many  engineers  fit  a  smooth  curve  to 
their  data.  Several  possible  physical  mechanisms  that  would  give  rise  to 
equilibrium  profiles  of  the  power-curve  type  have  been  described  (e.g.,  Bruun, 
1964;  Dean  1977).  Other  forms  that  sometimes  fit  profile  data  (e.g.,  log¬ 
arithmic,  parabolic,  etc.)  are,  like  the  power  curve,  everywhere  concave 
upward.  There  is  a  problem  inherent  in  the  use  of  such  curves  to  represent 
profile  response  to  changes  in  water  level;  e.g.,  Bruun  (1964)  found  that  the 
expression  for  offshore  deposition,  based  on  his  adopted  power  curve,  indi¬ 
cated  an  unrealistic  thickening  of  the  deposit  offshore.  Figure  23  and  the 
following  paragraph  show  why  this  and  similar  problems  occur  when  a  curve 
which  is  everywhere  concave  upward  is  adjusted  according  to  equation  (1). 


Figure  23.  Limitations  of  analytical  models  with  profiles  everywhere  concave  up.  Tangents 
will  have  only  one  point  of  intersection.  On  the  other  hand,  the  idea  of  ex¬ 
posing  a  trailing  edge  implies  that  the  offshore  slope  equals  Z/X,  i.e.,  the 
extension  of  the  offshore  slope  must  Intersect  the  profile  above  the  water 
surface.  This  is  impossible  if  the  profile  is  everywhere  concave  upward. 

Adjusting  an  equilibrium  curve  to  higher  water  leaves  a  trailing  edge  off¬ 
shore.  By  assumptions,  the  projection  of  this  surface  toward  the  shore  must 
intersect  the  profile  again  at  the  highest  point  of  wave  adjustment  (see  Fig. 
21),  but  the  tangent  of  any  concave -upward  profile  will  intersect  it  at  only 
one  point  and  everywhere  else  will  be  below  the  curve. 

If  it  is  assumed  that  the  offshore  slope  gradually  approaches  Z/X  near 
the  closure  depth,  then  concave  shapes  only  represent  the  inner  part  of  the 
active  profile.  Manipulating  such  curves  to  represent  adjustments  to  higher 
water  levels  inevitably  leads  to  unrealistic  consequences  offshore. 

d.  Inferring  Angle  of  Profile  Adjustment  from  Offshore  Slope.  As  dis¬ 
cussed  previously,  a  uniformly  sloped  trailing  edge  suggests  steady-state  con¬ 
ditions  (i.e.,  no  * signif icant  change  in  wave  climate,  profile  dimension,  or 
sediment  type).  In  such  cases,  direct  inference  from  slope  to  retreat 
(X/Z  ~  tan  a)  is  risky  because  forces  other  than  wave-induced  currents  may 
have  modified  bottom  slopes  over  the  long  timespan  of  profile  recession. 
Furthermore,  where  the  retreat  is  small  relative  to  total  width  of  the 
responding  profile,  the  mean  slope  over  this  short  distance  is  difficult  to 
measure  precisely.  Lastly,  errors  in  estimating  the  critical  depth  would  lead 
to  measuring  the  slope  at  the  wrong  place.  Nevertheless,  it  may  be  useful  to 
consider  the  types  of  geometry  implied  by  idealized  profile  adjustment,  com¬ 
pare  them  with  actual  profile  shapes,  examine  alternate  explanations  for 
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observed  shapes,  and  then  evaluate  the  results  in  light  of  all  ttie  other 
evidence  and  indications  at  hand. 


In  cases  where  the  only  data  are  hydrographic  surveys  (and  possibly  wave 
climate  but  no  data  on  backshore  deposits)  a  crude  first  guess  at  the  ratio  of 
retreat  to  submergence  could  be  made  directly  from  the  slope  of  an  apparent 
trailing  edge. 

Uniformity  of  slope  over  a  broad  section  of  the  critical  depth  on  adjacent 
profiles  is  striking  in  present  study  areas  (Fig.  24).  The  fact  that  this 
slope  equals  the  ratio  of  measured  retreat  to  lake  level  change  (x/z)  does 
not  necessarily  follow  from  the  observation  that  x/z  =  X/Z  (shown  in  Table  2), 
but  it  is  an  additional  observation  that  further  confirms  the  appropriateness 
of  the  sediment  balance  approach  for  long-term  predictions  on  the  Great  Lakes. 

The  regional  variation  of  profile  slopes  above  the  pinch-out  depth  (Fig. 
25)  reflects  active  processes  which  are  not  uniform  alongshore.  Rip  currents 
and  shoreline  undulations  are  other  expressions  of  such  lateral  variations. 
The  marked  divergence  of  slopes  below  the  critical  depth  reflects  deeper 
modern  processes  unrelated  to  surface  wave  action  or  relict  processes 
inherited  from  a  much  earlier  period  of  lake  evolution. 

V.  EXAMPLE  PROBLEMS 

The  following  problems  are  evaluated  on  the  basis  of  limited  amount  of 
available  survey  data.  They  provide  examples  of  the  basic  steps  in  applying 
the  proposed  method  for  profile  prediction.  If  these  predictions  were 
intended  to  support  actual  design  or  management  decisions,  a  more  careful 
evaluation  of  conditions  at  the  field  sites  would  be  required. 

***************  EXAMPLE  PROBLEM  1*************** 

G IVEN:  A  contemplated  change  in  the  regulation  plan  controlling  the  water 
supply  to  Lake  Ontario  would  raise  the  long-term  surface  elevation  0.3 
meter. 

F 1ND :  What  effect  would  the  higher  stages  have  at  the  eastern  end  of  Lake 
Ontario? 

ANALYS IS :  The  barrier  beaches  and  high  dunes  which  characterize  this  stretch 
of  shore  are  of  special  ecological  and  scenic  value.  Being  downwind  from 
the  major  storm  paths  across  Lake  Ontario  these  barrier  beaches  are  exposed 
to  the  highest  storm  waves  on  the  Great  Lakes,  but  because  of  relatively  low 
land  development,  few  protective  structures  exist  along  this  reach  of  the 
shore.  Sand  extends  lakeward  across  a  series  of  longshore  bars.  There  are 
no  known  rock  outcrops,  and  there  is  a  close  balance  between  southward  and 
northward  longshore  transport. 

EVALUATION  OF  TERMS: 


z  =  0.3  meter  Given 

h^  =  6.4  meters  From  site  16  (App.  G) 
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STATIONS  10-14 
1975 


STATIONS  15-19 
1975 


constanr  reference  slope  of  l  on  b9,  which  Is  based  on  the  regional  angle  of  effective  profile  ad justnent 
between  19b9  and  197  3.  The  maximum  and  minimum  estimates  of  the  pinch-out  depth  within  each  subregion  ar 
Similarity  among  bottom  slopes  is  evident  over  a  wide  area  near  the  pinch-out  dept  h,  — totu  i  nu<-.l 


13.4  me UTs 


Estimated  closure  depth  d  =  2.1  le  (  !  r.in  eij.  >) 

+  7.(>  meters  Average  IteigliL  ol  the  eroding  dunes  a t»o ve¬ 

st  i I lwate r  le ve  1  . 

=  21  meters  Sum  ot  the  two  values  obtained  above. 


Average  distance  ol  the  13.4-meter  depth  con¬ 
tour  t rom  shore.  The  vertical  ditui.i  should 
be  the  same  re  I  oroiiee  below  v:ii  rli  the  closure- 
depth  was  measured  in  the  previous  step. 

All  lit  the  material  eroded  Iron  the  upper 
beach  is  expected  to  remain  within  the 
bounds  ot  the  respond i  lg  proi ile. 

0.3(2,414)  l1 

— - - - - - =  34  meters  Evaluating  eiju.it  ion  (  1  ) 
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It  is  estimated  that  the  higher  stages  shift  the  equi  1  ibriau  snore  prat  il,- 
an  average  of  34  meters  inland  and  raise  it  0.3  meter  above  j>resent  coiidi- 
t  ions . 

***************  PXAMI’LK  I'KuliLIv!  J  **************  * 


Cl  IVES:  Assume  a  new  regulation  plan  is  proposed  to  noditv  the  in;  low  t  -  •  baru 
Michigan  and  Lake  Huron  via  the  St.  Marys  !<i  wr .  If  ldoptcd,  L  t .  i  -  pl.m 
would  lower  the  long-term  mean  surface  elev  it  foil  o'  hake  Michig.m  md  Luxe 
Huron  by  0.3  meter. 

KIND:  What  effect  will  the  lower  water  levels  have  on  shore  erosi-n)  -it  Liu. 

Indiana  Dunes  National  Seashore? 

ANALYSIS:  The  dredged  channel  and  navigation  structures  at  Micnigin  city, 

updrift  of  the  Indiana  Dunes  Nation,  il  Seashore,  block  some  >1  Lae-  potent  ia  l 
sediment  input  from  the  east .  Westward  longshore  transport  out  a  tue  dune 
area  thus  creates  a  sand  deficit  and  contributes  to  a  1  on.,-st  in  l  i  up  of. is  i.i-i 
problem  in  the  pa  rk .  As  lake  levels  tall  the  shoreline  withdriws  i.ni  tue 
beach  widens.  Assuming  take  currents  and  waves  ire  not  altered,  the-  :  ,.-.un 
to  reestablish  the  previous  prof i le  shape  at  a  lower  and  more  iukew-ird  posi¬ 
tion.  Longshore  losses  to  the  west  continue  to  exceed  the  net  supply  :r-n. 
the  east.  However,  offshore  where  the  bottom  slope  is  gradual,  lower!  p.  a 
the  water  surface  brings  bottom  sediments  into  a  shallower  tivdrnili.  rep,  i  me . 
This  results  in  ot t shore  sediments  moving  landward  to  steepen  nearshore 
slopes,  to  build  dunes  oil  the  widened  beach,  and  to  teed  t  n-  i  -  >n  c  -•  u- u  - 
currents  leaving  the  dune  area  to  the  west .  The  euiuul.it  iv.  .  I  I  u  t  s  .a  ii.es, 
.id  jus  tr, tents  can  be  estimated  using  equation  (I).  See  hands  -  i  'J  ■  M  !  1  a 
documentation  ot  shore  accretion  during  period  ot  do,  Uu..p  1  e„  l.v,  .. 

KVALLAT  ION  UK  I'KRMS: 

' .  i  veil 

average  I  rom  sites  2b  ind  24  (  \pp.  1  ! 

Dept  h  ol  protile  e  insure  =  1.1  n  tie,-,  .  > 

I. sli:nnled  averagi-  height  ot  dune-,  exp,-,  ted  : 

I  orm  oil  the  widened  beach  l.tkewirb  -C  tin 
present  l ore, tunc  . 


/.  =  -0.3  met e r 
it-  =  :>.)  me ters 

11.1  meters 
+  2.4  mi"  ters 


1  ->  meter 


Sun  hi  i no  two  valuer 


.ibt.i  i  loil  above  • 

X  =  3, 'Mu  meters  Average  distance  ot  Lho  1  1  .  1-meter  iimiuur 

J  r i jin  shore ,  basod  mi  lie  id  surveys. 

=  i  iM  1  s! tore  sands  are  expected  La  nave  onshore, 

and  L lie  wind  is  iu>L  expected  Li  carry  sand 
iiLand  past  tile  present  turedune. 

-  11M1 u> 


Lt  is  thus  estimated  that  Lowering  the  lake  Level  o.  3  meter  eiiecti.vl. 
shifts  the  aauil  !}i”iurrt  position  <>5  meters  Lakeward.  .Vs  discussed,  there  will 
still  he  a  net  loss  of  sand  due  to  net  transport  to  the  west;  there!  ore,  tne 
antuil  shoreline  Is  not  expected  to  advance  t>3  meters  lakeward.  A  reasonable 
interpretation  is  that  there  wi  l  L  be  a  long-term  gain  ot  to  meters  ol  hi  acn 
that  otherwise  wouLd  have  been  Lost  by  erosion  at  the  previous  water  levels. 
Dividing  b5  meters  by  the  appropriate  recession  rate  provides  an  estimate  o: 
when  the  avoided  erosion  would  otherwise  have  occurred. 

if  1 5  percent  of  the  offshore  sediments  are  in  the  clay-  ur  si  lt-si  «.• 

range  and  are  thought  to  he  too  fine  to  remain  in  the  active  shore  /.one, 

Lho  width  of  shore  saved  should  be  reduced  to  (1  -  0,13)  *  o')  meters  -  >  i 

meters.  Note  that  a  liberal  estimate  of  future  dune  heights,  L>,  would  also 
make  the  predicted  savings  more  conservative. 

ft*************************************** 

V  I .  SUMMARY 

Adjustments  ot  the  beach  and  the  nearshore  zone  to  long-term  change-,  in 
lake  level  were  monitored  between  1%7  and  N7o  along  a  5'1-kilometer  stretcu  .>: 
shore  centered  on  Little  Sable  Point,  Michigan.  The  bathymetry  of  this  region 
is  marked  by  a  sequence  ot  lour  to  live  longshore  bars  which  are  persistent 

from  year  to  year.  The  bars  are  continuous  over  tens  of  kilometers,  though  the 
pattern  is  disrupted  i.o  an  area  opposite  Little  Sable  Point.  In  the  cross  sec¬ 
tion  the  bars  are  much  less  regular  and  smooth  in  an  area  A  to  )  kilometers 

around  this  broad  protrusion  ot  the  shoreline.  The  bar-,  south  of  tile  point  are 
deeper  than  to  the  north.  drain  sizes  throughout  the  study  area,  both  on  the 
beach  and  along  the  bar  crests,  decrease  toward  tne  point.  The  longshore 
transport  converges  toward  Lhe  point.  The  tendency  ot  profiles  to  be  irregular 
in  areas  with  an  overabundance  ot  sand  lias  been  noted  elsewhere  bv  Uruun 
(l%d).  The  similarity  and  symmetry  ol  Lhe  other  patterns  suggest  a  common 

dependence  ot  all  the  discussed  variables  on  long-term  directional  wave 
ch  »rac ter  ist ics . 

The  longshore  drilt  which  converges  on  Little  Sable  Point  is  primarily  t rom 
areas  encompassed  by  the  study.  So  the  fact  that  the  longshore  transport  rates 
within  the  study  area  are  not  well  known  does  not  hinder  the  calculation  ol  a 
net  sediment  balance  tor  the  overall  region.  The  volume  ol  material  eroded 

from  the  upper  beach  over  A  years  averaged  A1  cubic  meters  per  me t e r- lengt h  ol 
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shore.  Zones  of  deposition  of  the  material  found  offshore  caused  the  cumula¬ 
tive  volume  measurement  made  along  the  profiles  to  increase  toward  zero.  The 
cumulative  volume  curve  typically  crossed  zero  and  displayed  high  positive 
values  before  settling  toward  a  constant  value  lar  from  shore.  There  was  no 
tendency  for  imbalances  on  one  profile  to  cancel  opposite  imbalances  on 
dire.ctly  adjacent  profiles.  However,  considering  the  overall  region,  onshore 
losses  closely  matched  offshore  gains;  there  was  no  evidence  of  significant 
exchange  beyond  the  surveyed  area. 

The  general  sequence  of  response  to  increased  water  levels  includes 
immediate  inundation,  gradual  migration  of  the  longshore  bar  sequence  up  the 
beach  slope,  and  increased  shore  recession  (but  at  a  rate  dependent  on  storm 
events).  Bar  migration  occurs  even  under  relatively  mild  wave  conditions. 
This  Lends  to  maintain  constant  bar  depths  even  while  the  mean  water  elevation 
is  changing.  Consequently,  the  barred  profile  becomes  compressed  toward 
shore.  The  erosion  of  shore  deposits  and  their  redistribution  lag  behind  the 
migration  of  the  bars.  Shore  recession  eventually  reestablishes  a  wider 
separation  between  inner  bars  and  the  waterline.  In  the  present  instance, 
reversal  of  the  lake  level  trend  occurred  before  all  the  material  deposited 
offshore  was  reshaped  to  reestablish  relief  on  the  outer  bar  comparable  to 
that  observed  at  the  beginning  of  the  study. 

VII.  CONCLUSION  AND  RKCOMMKNDAT IONS 

Hands  (  1979)  presented  a  set  of  shore  retreat  measurements  made  at  the 
present  profile  stations  over  various  periods  of  water  level  change.  The 
average  shore  retreat  for  a  given  change  in  water  level  was  approximately 
proportional  to  the  amount  which  t lie  water  level  had  risen  over  that  period  ot 
time.  It  was  suggested  that  this  linear  dependence  be  used  as  a  guide  in 
estimating  the  effects  of  future  lake  level  changes,  not  at  a  single  profile 
station  but  for  a  reasonable  stretch  of  similar  shoreline  responding  to  a  sim¬ 
ilar  submergence.  Qualitative  guidelines  suggested  how  the  estimates  should 
he  modified  to  reflect  differences  in  sediment  characteristics,  erosional 
forces,  and  the  length  of  time  considered.  Because  the  lake  level  and  shore¬ 
line  measurements  referred  to  changes  over  the  same  time  period,  no  allowance 
was  made  for  the  fact  that  the  shore  was  probably  out-of-phase  or  lagging 
behind  the  water  level  change.  It  was,  however,  pointed  out  that  some  lag  was 
inevitable  and  that  the  evidence  indicated  it  could  be  on  the  order  of  a  few 
years.  The  time  required  for  complete  readjustment  would  depend  on  the  energy 
available  for  sediment  redistribution. 

A  more  comprehensive  method  ot  estimating  profile  response  to  high  water 
is  developed  here  using  hydrographic  survey  data  to  extend  the  same  beach  pro¬ 
files  to  depths  of  more  than  12  meters.  A  simple  sediment  balance  equation 
predicts  the  amount  of  retreat  ultimately  necessary  to  reestablish  an  equilib¬ 
rium  profile.  Remackable  agreement  was  found  between  the  estimated  ultimate 
retreat  and  that  which  actually  accrued  3  years  alter  lake  levels  stabilized. 
Realistically,  the  equilibrium  model  also  overpredicted  shore  retreat  lor  the 
shorter  periods  of  sustained  lake  level  rise  before  stabilization. 

The  choice  of  whether  to  adopt  a  linear  relationship  between  retreat  and 
submergence,  making  the  qualitative  adjustments  as  discussed  in  Hands  (1979), 
or  to  apply  the  sediment  balance  approach  presented  here  will  depend  on  the 
tlmespan  of  interest,  the  amount  of  site-specific  data  available,  and  the 
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similarity  between  test  and  problem  conditions.  Initial  consideration  or  both 
estimates  may  be  useful  for  placing  high  and  low  bounds  on  the  expected  re¬ 
sponse.  Of  tile  two,  the  equilibrium  approach  is  more  objective  and  flexible 
as  it  takes  site-specific  character  ist ics  directly  into  account. 

both  approaches  must  be  qualified,  however,  for  the  possible  exposure  of 
nonsandy  substrates,  for  the  possibility  of  local  intervention  halting  erosion 
and  changing  longshore  balances,  and  for  the  conditions  that  existed  before 
the  period  of  application. 

Results  reported  here  are  promising.  However,  the  contrast  between  the 
extreme  simplicity  of  the  model  and  the  intractable  complexities  of  actual 
beach  and  nearshore  processes  emphasizes  the  need  for  careful  application  and 
further  evaluation  of  these  methods.  Careful  application  during  future  lake 
Level  cycles  should  provide  a  clear  indication  of  weaknesses  and  usefulness  of 
these  methods.  If  further  research  is  then  deemed  necessary,  the  section  of 
shore  studied  here  would  serve  as  a  good  test  site,  for  reasons  discussed  in 
the  text  and  because  of  the  available  past  record.  If  such  a  study  is  neces¬ 
sary,  plans  should  be  made  to  extend  the  study  over  at  least  a  full  cycle  of 
rising  and  falling  lake  levels  with  biennial  surveys.  The  full  longshore  unit 
should  be  covered  from  the  Ludington  pumped  storage  facility  to  the  White  Lake 
jetties.  The  positions  of  station  monuments  used  in  the  past  studies  are  we L L 
documented  and  can  be  reoccupied,  but  supplementary  profiles  should  be  estab¬ 
lished  between  these  situs  because  of  the  large  variation  in  shore  response 
observed  between  adjacent  stations. 

Because  of  the  required  profile  length,  boat  positioning  is  critical. 
Methods  other  than  the  usual  siting  on  surveyed  range  markers  are  necessary. 
The  time  of  year  for  surveying  is  also  important.  A  change  in  water  tem¬ 
perature  from  24°  to  13°  Celsius  within  2  hours  was  noted  during  sounding 
operations  in  this  study.  Extreme  temperature  changes  can  affect  the  repro¬ 
ducibility  of  soundings.  In  the  spring,  as  waters  warm,  a  sharp  temperature 
and  acoustic  gradient  develops  near  the  shore.  Significantly  cooler  water  is 
sometimes  trapped  in  a  series  of  pools  between  the  longshore  bars. 

Even  in  late  summer,  when  the  thermocline  is  typically  near  13  meters,  a 
change  in  wind  direction  can  quickly  flush  the  nearshore  zone  of  warm  water 
and  replenish  a  series  of  Longshore  pools  with  cold  bottom  water.  Calibration 
of  the  sounding  instrument  with  a  reflector  suspended  over  the  side  of  the 
boat  (a  bar-check)  should  be  done  from  the  surface  to  ttie  maximum  profile 
depth,  and  in  a  water  column  essentially  like  that  at  the  profile  site. 
Variations  which  cannot  be  eliminated  by  field  adjustments  can  be  corrected 
during  data  reduction  stages  if  careful  notes  are  kept  of  bar-check  results. 

A  good  evaluation  of  the  simple  profile  response  model  presented  here 
should  be  relatively  easy  after  the  lakes  have  undergone  another  long-term 
cycle.  However,  obtaining  the  field  data  to  significantly  improve  the  situa¬ 
tion  would  be  a  more  difficult  and  expensive  undertaking. 
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APPI'NUH v  H 


rM.D'IAT  I'l'is 

This  appendix  presents  net  nrnftle  changes,  cumulative  volume  curves 
n’Vn,  and  average  thicknesses  of  these  volumes  ( Th  1 .  The  calculation 
procedures  are  described  in  Section  Til,  h,  h  and  illustrated  in  Figure  l"1. 
Table  1  in  the  text  tabulates  the  rates  of  chance  ip  cubic  meters  ner  meter 
ner  vear. 
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APPENDIX  C 


VARIATIONS  IN  EXPOSURE  TO  STORM  WAVES  AT  THF  ORFAT  I  .AXES 

This  appendix  presents  location  maps  (Figs.  0-1  to  C-S)  and  wave  height 
measurements  (Tables  C-I  to  C-5)  for  the  Oreat  Lakes  study  sites.  Heights 
along  the  fence  shown  below  are  related  to  the  energy  of  storm  waves  obtained 
from  a  numeric  hlndcast  model  hv  Resto  and  Vincent  (197ha,  197hb,  197f>c,  1977, 


Table  C-l.  Wave  heights  fin  meters)  at  Mke  Superior  sturlv  sites 
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47.11 

88 . 98 

5.4 

29 

Rl  PR  1  01  !  ,  Ml 

47.24 

88 . 77 

5 . 3 

30 

CALUMET,  Ml 

4  7. 30 

88 . 55 

i 

31 

1  Alil.I  p  1  VI  1 1  ,  Ml 

47.53 

60 . 34 

4 . 9 

32 

EAGLE  HA  k  NOR .  Ml 

47.53 

88.14 

4  .  7 

33 

COPPf  R  IIAkl  llk,  Ml 

47.53 

8  7.91 

4.5 

34 

SCHLATTER  LAKE,  Ml 

4  7.53 

87 . 70 

4 . 3 

3S 

MON  1  T  III  1  ISLAND,  Ml 

47.38 

67.49 

5.4 

36 

KtUEF  NAU  PI  1 1  N  I  ,  MI 

47. 3H 

87 . 70 

c  « • 
>J  .  vj 

37 

PUINJ  ISAPII.1.1,  Ml 

4  7.24 

87.91 

5.4 

38 

TRAVERSE  Pi  1 1  N  I  .  Ml 

47. 09 

88.14 

5.5 

39 

PIOIIAMINI.,  Ml 

46 . 95 

08 . 35 

4 . 7 

40 

HURON  kIVfk  POINT, 

MI 

46 . 95 

87.91 

4 . 4 

41 

MU  PAY,  Ml 

46.95 

87.70 

4.8 

42 

GARLIC  POINT  ,  Ml 

46.81 

87.50 

4 . 9 

43 

MARiJUl  Til,  Ml 

46.66 

87.29 

5 . 1 

44 

PIT  R  ION,  Ml 

46.66 

87 . 08 

5.2 

45 

AU  TRAIN  PAY,  MI 

46.66 

86 . 86 

5.1 

46 

GRANO  ISLAND,  Ml 

46.66 

86 . 65 

5.3 

47 

GRANJi  POM  ML  POINT, 

Ml 

46.65 

86.45 

5.4 

48 

AIJ  SAPLE  POINT,  MI 

46 . 79 

86.23 

f  f 

>J  •  J 

49 

GRAND  MARAIS,  MI 

46.78 

86.02 

5  .  Ti 

50 

SUCKER  RtVCR,  MI 

46.78 

85 . 80 

5 . 5 

51 

DOR  PARK,  MI 

46 . 70 

85.60 

5.6 

52 

LITTLE  LAKE  HARPOR, 

MI 

46.77 

85 . 39 

5.6 

53 

CRISP  POINT ,  MI 

46.91 

85.18 

5.3 

54 

PARADISE,  Ml 

46 . 62 

84.97 

3.0 

55 

POINT  IROIJOniS,  MI 

46.46 

84 . 75 

3.9 

56 

GROS  CAP,  MI 

46.61 

84 . 76 

3.0 

57 

GOULAIS  RIVfR,  MI 

46.76 

84 . 75 

5.5 
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Table  0-2 

.  Wave  heights  (In 

meters ) 

at  Lake 

Michigan 

study  sti 

SITE 

UK  AT  I  no 

l AT l TOOL 

LONG  I  tOIlF. 

h5 

I 

Mt'EM  tin  !  A*  .  hT 

45 . 64 

85 . 37 

5.4 

n 

HAUl.HK  SPRING,  nt 

45.50 

85 . 32 

3.7 

3 

1  1SMI  KHAN  I'.LAno,  Hi 

45.3  7 

85.50 

4.3 

4 

F  RAVI  USE  I'M.  Hi 

45 . 23 

85 . 72 

4  .  1 

*» 

TRAVI  I'Ll  till,  Hi 

45.10 

86.15 

3.7 

6 

HAN  llliU.  Hi 

44.95 

86.15 

3.9 

? 

HI.  A 1  T  t  If.sE,  HI 

44  .  ao 

86.17 

4.6 

a 

FRANRHlKI.  HI 

44 . 60 

86  .  W 

5.2 

9 

ARCAG 1  A  ,  HI 

44.53 

86.40 

5.2 

10 

(INF  K  AH  A,  Hi 

44 . 30 

8  ,  40 

4 . 3 

1  1 

MAN151H  ,  Hi 

44 . 24 

06.41 

4 .5 

12 

HUi  GAI-LE  HIIINF.  HI 

44 . 09 

86.63 

4  .  I 

13 

LUPINUIHN,  Ml 

43.94 

86 . 64 

.  3 

14 

Hf.HIUATFK.  HI 

43 . 80 

86.66 

5. 1 

ir. 

i.nn.i  saiu  pniNi 

HI 

43.65 

86 . 66 

5 . 4 

16 

HF  NONA  .  H  1 

43.52 

86 . 67 

5.4 

17 

HUNT  Alii  It  ,  nl 

43.36 

86 . 4R 

5.7 

18 

MHSKF  flHN  ,  Hi 

43.23 

86.50 

5.7 

19 

CRANK  HA VI  N  ,  HI 

43.06 

86.32 

5.6 

20 

ISRANUF'.  RAHlIiS,  Hi 

42.93 

86.33 

5.6 

21 

HULl.ANli,  HI 

42. 78 

86.33 

5.3 

22 

tKllK.t.Ar,,  HI 

42.64 

86.35 

5.4 

23 

SOll  1 H  MAV1N,  Hi 

42.48 

86 . 36 

5 . 3 

24 

S  SOUTH  HAVEN,  HI 

42. 34 

86.37 

5.2 

25 

HINTON  HAKLUR,  HI 

42.21 

86 .57 

5 . 0 

26 

S  ST,  JOSEPH,  HI 

42.06 

86 . 58 

5 . 0 

27 

Nt  U  MU  1  ALO.  HI 

41.93 

86. 7H 

4.8 

28 

MICHIGAN  CITY,  1  NO 

41  .'78 

86.99 

5.2 

29 

RJRNS  HARHOR,  1NU 

41.79 

87. 18 

5.4 

30 

CHICAGO,  11.1. 

41.80 

6  7.38 

5. 1 

31 

CHICAGO  SHIP  (ANAL., 

IU 

41.95 

87.56 

5'  .  4 

32 

EVANSTON,  ILL 

42.10 

8  7 . 56 

5 . 3 

33 

HIUHl.ANl'  f'AKK,  ILL 

42.26 

87.73 

5.2 

34 

UAHKFGAN,  ILL 

42.40 

8  7. 73 

3.8 

35 

KENOSHA,  UT 

42.54 

87.73 

5 . 2 

36 

S  RACINE,  UI 

42.69 

87.71 

3.9 

37 

N  RACINE,  UI 

42.83 

87.70 

5.2 

38 

S  HILUAUKrC,  UI 

42 . 97 

8  7.69 

5.0 

39 

milua(ikf:e,  ui 

43.12 

87.68 

4.8 

40 

S  PORT  UASHTNbfON,  UI 

43.27 

87.68 

3.9 

41 

PORT  U ASH I NO T ON ,  UJ 

43  41 

87.67 

3.7 

42. 

N  PORT  WASHINGTON.  UI 

43.55 

8  7.66 

3.8 

43 

S  SHF  80YGAN,  U T 

43.69 

87.65 

4.0 

44 

N  SHE60Y0AN,  UI 

43.84 

87.65 

4 .7 

45 

MANITCIUUC,  UI 

43.98 

87.64 

3.8 

46 

TUO  RIVERS ,  UI 

44 . 13 

8  7.43 

3.7 

47 

RAUt.F  Y  KJJNT,  UI 

44.27 

87.42 

4 .5 

48 

KF.  UAIINf’E  ,  UI 

44,42 

87.41 

3.8 

49 

Af.  cum  A,  UI 

44.56 

87.20 

3.6 

50 

STURGEON  HAY  CANAL, 

UI 

44.70 

B7 .20 

3.6 

51 

JACKSONPtlRl  ,  UI 

44 . 84 

87. 18 

4.7 

52 

RAIIEY3  HAKPtjR,  UI 

44.98 

86.98 

4  .5 

53 

N  CANAL  1.U1H1,  UI 

45.  J4 

86.96 

3.6 

54 

WASHINGTON  JVLANt', 

UI 

45.27 

86.74 

3.4 

55 

E ISHI WHAN  SHOAL.  UI 

45.41 

86 . 73 

3.2 

36 

ESC ANAPA ,  MJ 

45.54 

86  52 

4 . 7 

37 

POINT  AUX  EARiJMES, 

HI 

45.68 

86 . 30 

4 . 1 

58 

N  POINT  AUX  PARIJHES 

,  HI 

45.83 

86,29 

3.3 

39 

HANTS T 1 01 II  ,  HI 

45.83 

86.08 

3.9 

60 

PORT  ISLANIi,  HI 

45  81 

85 . 88 

5.7 

61 

POINT  PAIIIRSON,  HI 

45.94 

85 . 66 

2.6 

62 

H fl.i.E  rf'UHFNS  RIFE, 

HI 

45.94 

85.45 

n  **> 

63 

SAIILF  ST F  HAR 1 1  ,  HT 

45.93 

85.25 

2.4 

64 

wreuoor r  lakf:,  mi 

45.92 

85,04 

5.7 
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Table  C-3.  Wave  heights  (in  meters)  at  Lake  Huron  study  sites. 


SHE  LOCATION 


LATITUDE  LONGITUDE  h 

5 


J. 

PURI  HURON,  MI 

43.0  2 

82 . 33 

LAKE PORT,  MI 

43. 16 

82 . 33 

3 

LEXINGTON,  MI 

43 . 31 

82 . 33 

4 

PORI  SANILAC,  MI 

43 . 45 

82 . 52 

C 

w.l 

FOREST VI LI  ,E  ,  MI 

43.59 

82 . 52 

6 

HELENA ,  MI 

43.74 

82 . 52 

•*v 

( 

HARBOR  BEACH,  Ml 

43.88 

82 . 52 

8 

HURON  CITY,  MI 

44 . 0  3 

82.71 

9 

PORT  CRESCENT,  MI 

44.03 

83 . 1 1 

10 

ENTRANCE  SAGINAW  BAY,  MI 

44 . 0  3 

83 . 30 

1  1 

I  A WAS  CITY,  MI 

44 . 1.7 

83 . 30 

1 2 

OSCODA  All  SABLE,  Ml 

44 . 31 

83 . 1 1 

1.3 

GREENBOSH,  MI 

44 . 46 

83.11 

14 

HARR ISV.IL I..E  ,  MI 

44 . 60 

83.11 

15 

BLACK  RIVER,  MI 

44 . 74 

83 . 1.1 

16 

OSS I NEK E,  MI 

44 . 89 

83.11 

17 

NORTH  POINT,  MI 

45.04 

83 . 1 3. 

18 

ROCK  PORT,  MI 

45. 18 

83.30 

1.9 

STONE  PORT,  MI 

45 . 32 

83.30 

20 

ADAMS  POINT,  MI 

45.46 

83.50 

21 

ROGERS  CITY,  MI 

45 , 46 

83.70 

/;> 

HAMMOND  BAY,  MI 

45 . 61 

83 . 90 

23 

CORD WOOD  POINT,  MI 

45 . 75 

84 . 09 

24 

POINT  DOLOMITE,  MI 

45 . 89 

84 . 09 

25 

DEI  OUR  RE EE ,  MI 

45 . 89 

83 , 90 

26 

W  END  DRUMMOND  ISLAND,  MI 

45.75 

83 . 7  0 

'"i  *7 
*•:.  < 

FALSE  DETOUR  CHANNEL,  MI 

45.75 

83 . 50 

28 

COCK BURN  ISLAND,  MI 

45 . 75 

83.30 

6 . 2 

5 . 8 

4 . 7 
5 . 0 

3 . 7 

5 . 4 

5 . 7 

5.8 

5.9 
6. 0 
6  . 


6,2 
6.1 
6, 1 
6,1 
6 , 0 

5.2 

4 . 6 
5 . 0 

4 . 7 

4 . 2 

5.4 
5 . 0 

4 , 9 

5 . 4 
5 . 6 


i 


NO  -o  m 


Figure  C-S.  Study  sites  at  Lake  Ontario* 


Table  05.  Wave  heights  (In  meters)  at  Lake  Ontario  study  sites 


LOCATION 


LATITUDE  LONG  I  TU  0,.:: 


FORT  NIAGARA ,  NY 

43 . 4  3 

79 , 0  3 

WILSON ,  NY 

43 . 43 

78 . 83 

APPLETON ,  NY 

43 . 4 'j 

78 , 63 

THIRTY  MILE  POINT ,  NY 

43 , 45 

78 , 4/ 

LAKESIDE  PARK,  NY 

43 . 45 

78 . 25 

KENDALL,  NY 

43 , 47 

78 , 0  " 

NORTH  HAMLIN,  NY 

43 , 47 

77 , 8 7 

WEST  ROCHESTER,  NY 

43 . 48 

77 . 65 

EAST  ROCHESTER,  NY 

43.37 

77 , 47 

POL TNEYVILLE ,  NY 

43 . 37 

77.23 

SODUS,  NY 

43 . 38 

77 . 07 

WOLCOTT,  NY 

43 . 38 

76 . 87 

FAIRHAVEN,  NY 

43 , 48 

76 . 61 

OSWEGO,  NY 

43 , 63 

76 . 5) 

LAC] NA ,  NY 

43 , 63 

76 . 23 

BELLEVILLE,  NY 

43 , 77 

76 . 23 

GALLOO  ISLAND,  NY 

43 .91 

76 . 5) 

